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ABSTRACT 


This thesis considers the problem of determining 
optimal product mix and machine loading for single and multi- 
stage production systems. It is assumed that the Group 
Technology (GT) concepts have been utilized to classify the 
jobs under consideration for manufacture on the production 
systems into groups. Each job has a lot size comprising of 
identical units and each production stage has only one 
machine. In case of multi-stage production systems* it is 
further assumed that the job operations do not have any 
technological ordering. Corresponding to each stage, the 
job available for processing either have prespecified 
production times or the production times are to be determined 
optimally considering machining speeds as decision variables. 
The optimization is to be performed considering maximization 
of production rate or maximization of total profit. 

For each one of the above stated cases, an optimizing 
algorithm which uses the concepts of branch and bound proce- 
dure is developed. 

In order to save on the computational effort involved 
in the use of branch and bound procedure, for each case a 
heuristic method has been developed. The heuristic method 
helps in reducing the size of the problem to be handled by 
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the branch and bound procedure. Illustrative examples are 
given to explain the various steps of the optimizing and 
heuristic algorithms. 

The optimizing and heuristic algorithms have been 
computerized and implemented on DEC-1090 Computer System. 
Based on the limited experience of solving the numerical 
examples? it was observed that the heuristic algorithms have 
computational superiority over the optimizing algorithms. 



CHAPTER I 


INTRODUCTION 


1.1 GROUP TECHNOLOGY : 

Group Technology is a very progressive method of 
organising production and especially it is becoming very 
popular in those industries which are engaged in medium and 
small batch productions. Group Technology concept identifies 
and makes use of sameness or identical nature of operation 
processes and parts in design and manufacture. Similar parts ' 
are grouped together into one part family based on similarity 
of their design and dimensions, geometrical shapes and technolo- 
gical requirements, and according to the needs to process 
designated part families, machine group or cells are formed. 

1.2 REVIEW OF PREVIOUS WORK! 

Analytical methods suggested by Burbyidge [l] take into 
consideration those informations which are contained in route 
cards. However, he does not consider design and geometrical 
similarit 3 >-. Batra and Rajgopalan [2] also use informations 
given in route card for each job and using graph theoretic 
approach they form machine cells. McAuley [3] uses similarity 
coefficient for every pair of machines to form group of machines. 
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4 

El-Essawy[4] does component flow analysis to form machine and 
component grouping. 

Many researchers have proposed models to schedule 
products on machines available in machine cells. Petrov [5] 
has developed four interrelated scheduling models for different 
type of route sequences and component flows. This approach does 
not guarantee optimality mathematically. Ham and Hitomi [6] have 
suggested a technique to schedule various grouped products on 
multi-stage production system. 

In addition to the problem of group formation and group 
scheduling, machine loading and product mix decisions represent 
major problem areas for group production planning and schedul- 
ing. Machine loading problem determines the total number of 
units of various jobs to be processed within limited capacity 
of various production facilities available, while product mix 
decision determines the kinds of jots to be processed within 
limited production time. Hitomi and Ham 1.7, 8] have considered 
the problem from the view point of Group Technology for a single 
machine production s^^stem. They have presented models for' two 

distinct situations. In the first case, the unit production 

i 

time is prespecified [7] 1 while in the other case, it is 

C 

dependent on the machining speed which is treated as a decision 
variable [8], They maximize total number of units produced of 
various products and in process do not account for weightages 



or importance of various products to be produced. Moreover, the 
experience of the author in using Hitomi and Ham' s algorithm on 
a number of problem suggests that their algorithm does not 
yield optimal solution necessarily as claimed b 3 r them. 

1.3 SCOPE OF PRESENT STUDY s 

Basically, this thesis is concerned with the machine 
loading and product-mix decision problem considered by Hitomi 
and Ham [7,8] with added features to make the problem more 
realistic. In real life,single stage production environment is 
not very common. Further, besides maximizing the total amount 
of production, the management may desire to evolve a product mix 
and machine loading strategy which maximizes the profit. There- 
fore, in the present work, an attempt is made to develop models 
and solution methodologies for multi-stage situation with 
maximization of total production or maximization of profit as 
the criterion of optimization, considering the unit production 
time as prespecified or dependent on the machining speed ( a _ 
decision variable ). In general, the solution methodologj? 
exploits the branch and bound concepts for its development. 

Since, the use of branch and bound approach involves con side- ; 
rable effort, heuristic based computationally efficient solution j 
methodologies which initially reduce the size of the problem 

I 

and then find the solution applying branch and bound concept j 
to the largely reduced problem, are also presented. ; 
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Further, the algorithm for single machine loading 
problem and product mix decision of Hitomi and Ham [7,8] has 
been modified. The modified algorithm ensures optimal solution 
to the problem. 

1.4 ORG AIII SAT I ON OF THESIS: 

The cases of single stage production system have been 
dealt in Chapter II. In this chapter besides the modified 
version of algorithm of Hitomi and Ham [7,8], optimal and 
heuristic solution procedures for maximization of profit are 
also presented. Various variations of the multi-stage produc- 
tion system have been considered in Chapter III. Numerical 
examples are presented to illustrate the solution methodologies 
proposed for each of the variations of the single and multi- 
stage production systems considered. In Chapter IV, conclusions 
based on the present study alongwith suggestions for further 
work are presented. 



CHAPTER II 


SINGLE STAGE PROBLEM 

2.1 STATEMENT OF THE PROBLEM s 

Consider a single stage production system with limited 

available production time.N types of jobs are available for 

production. Based on part-family concept of group technology, 

these jobs are classified into N groups. Let groups 

6 

G^ (i = 1,2,..., N ) contains Ih jobs. 

The objective is to determine an optimal machine loading 
and product mix such that within the available production time, 
the quantities produced of each job do not exceed the lot size 
limit. The optimization is carried out to maximize either the 
production rate or the total profit. 

Corresponding to each group, there is a group production ; 

t 

time which is comprised of the group set-up time and the sum 
of the job production times of the jobs in this group. The job 
production time includes the job set-up time and the time 

required to produce the various units of that job. The unit : 

I 

production time required to produce single unit of any job may 
be prespecified or dependent upon the cutting speed (a decision | 
variable) . 
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2.2 ASSUMPTIONS? 

The various assumptions are; 

i) The single stage production system contains only one 
machine . 

ii) Only one unit of the job can he processed on the machine 
at a time. 

iii) A limited amount of time is available for production on 
the machine. This is the actual productive time available 
on the machine and as such do not include time lost due to 
breakdowns, etc. 

iv) Each unit of the job must be processed to completion. 

v) All the units of a job are produced together. Similarly, 
all the jobs belonging to same group are produced together. 

2.3 NOMENCLATURE % 

Following notations are used for the development of the 
mathematical models: 

a. . preparation time for job J. . (min/pc) 

-i- J 1 J 

bfj tool replacement time for job (min/pc) 

B. . profit per unit of job J. . (Rs./pc) 

C. . l-minute tool life machining speed for job J. . (m/min) 

-^d i 

C u total capacity of the machine utilized (min. ) 
d available production time (min) 

E. . high-efficiency machining speed range for job J. . (m/min) 

d IL ."1 i 
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G i 


N ± > 


J id 


10 

n id 


N 

N. 


N. 


P 


i0 


(c) 


i-th group or part family 
group index ( = 1,2,..., N a ) 
job index (= 1,2,.. 
j-th job in i-th group 
lot size of J. . (pcs) 

slope constant of Taylor tool life equation 

total number of jobs available for production 

total number of part-families or groups 

total number of jobs in group G^ 

unit production time for J. . (min/pc) 

i<D 


p. . minimum-production -cost unit production time for job 


il 


(t) 


(min/pc) 


p. . ' maximum-production-rate unit production time for job 


il 


J. . (min/pc) 
X J 


p. . (v. .) cutting speed dependent unit production time (min/pc) 
xj x J 

P. . job production time for J. . (min) 

i j 10 

P- X maximum-product ion-rate job production time for J. . (min) 

1 j X 

q. .(v. .) speed dependent unit production cost (Rs/pc) 

XJ X J 


S id 


10 

T. . 
10 


U 


V id 


job set-up time for J. . (min) 

X J 

Group set-up time for group G^ (min) 
actual machining time for J. . (min/pc) 

1 J 

tool life for J. . (min/edge) 

X J 

set of all parts available for production 

1,2,..., N } 


- h ; i 


{ J ±d / 0 — 1,2,... 

machining speed for J . . (m/min) 

X J 
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v . . 

V i0 

V i0 

w. . 
10 

X i0 

X i 


Y 

Y' 

z 

z 

z 

z. 

a 


(c) 

(c) 

(t) 


minimum-production -cost machining speed for J. - 
high efficiency machining speed for J. - (m/min) 

-L J 

maximum-production rate machining speed for J. 


10 


max 


^iO 

Y . . 
10 

6 

Aid 

0 

n 

* 


quantity of J. . to be produced (pcs) 

J 

0-1 type variable for J. - 
J 10 

i 

0-1 type variable for group 

total production cost (Rs« ) 

total variable production cost (Rs.) 

value of objection- function 

lower bound on objective function value 

upper bound on objective function value 

maximum value of objective function 

direct labour cost and overhead (Rs./min) 

machining overhead for J . - (Rso/min) 

10 

tool cost for J. . (Rs./edge) 

-L J 

slack time (min) 
machining constant for J . . 

-L J 

null set 

lagrange multiplier 
shows optimal value 


Following additional notations will be used in 


(m/min) 

(m/mm) 


solution procedure i 
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D consumed production time (min) 

n all nodes belonging to N1 

N1 set of all descendants 

R set of remaining jobs at any node 

s time still available for proiuction (min) 

S set of selected jobs at any node 

shows candidate for acceptance or rejection 
[A] is a Guassian notation imploring a greatest integer less 
than or equal to A. 

2.4 MATHEMATICAL MODELS AND OPTIMAL SOLUTION PROCEDURES* 

Mathematical models for the various cases of the single 
stage problem stated in Section 2.1 are presented alongwith the 
optimal solution methodologies. 

2.4.1 Ca se 1; Maximization of Total Units P roduc ed* 

In this case we maximize the total amount of production 
irrespective of the relative values associated with the various 
jobs to be processed. Two variations of the basic problem are 
considered. In the first variation we assume that the unit 
production time is prespecified while in the other case we 
assume that the production times are dependent upon the machining 
speeds. The machining speeds are also to be determined opti- 
mally considering the workpiece - cutting tool-machine tool 
considerations. 



2.4. 1.1 Model Is Unit Production Time is Pre specified s 
Constraints ; 

Various constraints on the problem are listed below, 
i) Constraint on Machine Utilization; 


This constraint ensures that total time consumed for 
production does not exceed the available capacity of the machine, 
Mathematically, this is represented ass 


N g 


1 


E (S i X i * * *ij> 

= 1 3 - L 


( 2 . 1 ) 


where . 


P^j is job production time for producing w^ units of 


job Jb ^ . It is expressed as, 


p id 


s. . + w. . p. . 
10 10 *10 


( 2 . 2 ) 


_ (1 if job Jj_j is selected for productions/- ,\ 

lj “ ‘0 if job J. . is not produced 

J 


and, 


x i 



if 


1 

£ x. . > 1 
j=l “ 


otherwise 


(2.4) 


ii) Constraint on Units Produced: 

The quantity w. . to be produced of job J. . should not 

lj 1 

exceed the lot size specified for it. Thus, 
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O bjective function : 

Since the objective is to maximize the total number of 


the units produced, we maximize a function presented as; 

N N. 
g i 

Z = I I w. . x. . 
i=l j=l 


( 2 . 6 ) 


Solution procedure ; 

Let Q represent the time needed to produce all the jobs 
to their lot sizes. Then, 

N N. 

1 l 

(2.7) 


Q = (Si + ^ (s.. + 1,, P,,)} 


-‘-J il 12' 

If Q < d,then all the jobs will be produced equal to their lot 
sizes and the optimal solution will be: 


where. 


■¥r 

X i 2 

= i 

¥ J. . e U 

12 

(2.8) 

w? . 

= i. . 

12 

¥ J ij £ U 

(2.9) 

Z* 

= Z 

max 


(2.10) 

z 

max 

N N. 

g i 

= I E 

i=l j=l 


(2.11) 


If Q > d, then some of the jobs are either not accepted for 
production at all or are not produced to their lot sizes. In 
such a situation, obviously, 


Z* < 


max 



12 


A branch and bound procedure is developed to optimally 
ascertain the 'machine loading and product mix for such a 
situation. 

In the branch and bound procedure, at any stage we branch 

a node n for which the upper bound value Z, is maximum among 

all descendants. Z is calculated as the sum of the lower 

S, (1-X. )+s., 

bound value Z and maximum of [ s/ { — — ^ J- + p. . }] among 

id 13 

all remaining dobs at that node? where s shows the remaining 
time available for production. Initial^, Z = 0. If at this 
branching d°h Jhf-' is selected for production, the lower bound 
value Z, is incremented by wrr. 

Thus, the branching procedure finds the lower bound and 
the upper bound values for each of the descendants. The node 
with maximum upper bound value is identified. For this node, 
if the upper bound value is equal to the lower bound value, 
then branching procedure is terminated and this node carrying 
the maximum upper bound value gives the optimal solution. 

Optimizi ng algori thm s 

The following is the step by step procedure for deter- 
mining the optimal machine loading and product mix. 

Step I s Calculate Q, from (2.7). If Q < d, then go to Step 2J 
■otherwise, go to Step 3. 

Step 2 ; Manufacture all the jobs to their lot sizes. The optimal 
solution is given by the following relationships. 
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<3 - w « - h 3 ¥ J 13 £ u md Z * ■ W 
Terminate the procedure. 

Ste p 3} Follow Branch and Bound Procedure. ; 

Select an initial node; set n = 1, N1 = {n}t For this 

node n, set R = U, S = 0 9 s = d, Z = 0 and = 1 ¥ i. 

Find a nob Jr? from the nobs J. . s R for which 
u il 

(St + Stt)/1t-? + p-r-r is minimum. Find 
i io y/ ii ^ii 

Z = [d/{(Sj + + Pjj}] a nd go to Step 4. 

Step 4 ; Find a node from the nodes n e N1 for which Z is maximum. 

Denote this node by n and J4-r in n by Jrr. For the node 

J il J il 

n, if Z s Z, go to Step 8 j otherwise, go to Step 5. 

Step 5 : Find D = D + S~X~ + src- + wc"r p~~, where wo"v = max w„. A $ 

11 11 11 * 1 ! il il 9 

w. . is varied from zero to lr~ till D <_ d. If wrv = 0, 
go to Step 65 otherwise, go to Step 7. 

<v 

Step 6 ; Remove the job Jr-v from the set R of the node n. Set 

-*-U 

R = R~{ Jc'o'} . If R = 0, set Z = Zj otherwise, find 
S.X. + s. . 

min ( — — + p ) and denote this Job by JVr. 
J i;j eR ii 10 13 

Set Z = Z + [s/ftSjXj + Sj-)/lj^ + Pl -B. 

Return to Step 4. 

, j t\j 

Step 7 s Branch the node n into two nodess 

First nodes Set n = n+1 and N1 = N1 + {n} - {n}. Then 
for the newly created node n, R = R - { Jrq'} , S = S, 

D = D, s = s, and Z = Z. If R = 0, set Z = Z; otherwise, 
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find min {(S.X. + s. .)/!.. + p. .} and denote the 
J. .eR 11 ID iJ 10 
10 


corresponding job by Set ^ = Z + + sj-r)/l-, 

1 J X 


+ P I3 }] ‘ 


Second node: Set n = n+1 and N1 = N1 + {n}. Then for 
the neyly created node n, R = R, S = S + { Jon} , D = D, 

X J 

s = d - D, Xo' = 0 and Z = Z + won. If R = 0, set 

i ““ x o 

Z = Zj otherwise, find min {(S.X. + s . . )/l. .+p . - } . 

t P p XX X J 1 j X J 

id 

Denote the corresponding j°To by and set 
z = z + [s/USjXj ♦ SI3 )/1 I3 + P I3 }]. 

Return to Step 4. 


St eps : Optimal solution will be contained in node n, Set, 


,* 




3 = S, Z = Z and stop. 


The value of won' required in Step 5 is determined 

X J 

using the following relationship. 


w'y'Y 
• 10 


d-D - SrXr 
xx 

prcr 

10 


- sc 


iX] 


If won > Ion, we set won = Ion. 
10 10 10 10 


Nu merical exampl e; 


The example given is based on the data of a single 
stage problem considered by Hitomi and Ham [?]. 

Consider a single stage problem where ten jobs are 
available for processing. These jobs are divided into four 
groups. The relevant data for the problem are given in Table 2. 
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Table 

2.1s Basic 

; Data for 

Numerical 

Example 1. 

GROUP NO. 

PART 

LOT SIZE 

GROUP SETUP 

JOB SETUP 

UNIT PRODUCTION 


NO. 

pcs 

TIME (min) 

TIME 

(min/ lot) 

TIME (mm/pc) 


1 

30 


19 

6 

3 

2 

40 

40 

8 

2 


1 

20 


10 

17 

2 

2 

50 

35 

9 

9 


3 

30 


15 

7 


1 

20 


5 

12 

3 

2 

60 

20 

13 

6 


1 

20 


6 

16 

4 

2 

10 

45 

10 

15 


3 

40 


20 

13 


Available Time d = 40 hours. 


Solution s Using the various steps of the algorithms, we 
obtain, 


lz Q = 5105 mins. 

Since Q > d, go to Step 3. 

Step 3 * n = 1, N1 = {n} = {l}. For node 1, set 

R = {J 11 » J 12 * J 2 i» J 22 J J 23’ J 31 * J 32’ J 4l* J 42’ J 43^ 
S = 0 , a = 2400.00 mins. , D = 0, Z = 0 and = 1 ¥ i. 

Calculate {(S. •*- s. . )/l. . + p. .} for J. . e R. The 

-L nj 1 «j -*-vJ 

minimum value of 6.55 min/pc corresponds to job J^ 2 * 
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Step 4 : 




Thus, Z = 366. Proceed to Step 4. 

Since for node 1, Z ^ Z, go to Step 5. 

For node 1, D = 393 mins, and w^ = 60. Therefore, 
go to Step 7. 

Branch out from node 1, two nodes, 2 and 3- 
At node 2: Set n = 2, N1 = {l} + {2} - {l} = {2}. 

For the newly created node 2, set R = R - {J^ 2 } 

= (J-q* J 12 » J 2i’ J 22 9 J 23 9 J 31 9 J 41* J 42* J 43^ 

s = 2400 mins., D = 0 and S = 0. Since R £ 0, 

calculate (S^ + s..)/-!^ + P ±3 for J.. a R. 

Identify the job for which the value of this expres- 

sion is minimum. For the present problem the minimum 

value of the expression corresponds to job J-^. 

Similarly, at node 3: Set n = 3, N1 = {2} + {3} = {2,3 

For the newly created node 3, set 

r = r = (J-Qy J 12 > J 21 + J 22 9 J 23 9 J 31 9 J 41 sJ 42 ,J 43^ 
D = 393 mins., s = 2007 mins, S = { J^ 2 1 > Z = 60 and 

= 0. Since R £ 0, we find the minimum of 

(S-X. + s. .)/!.. + p . . . Its value is 6.8 min/pc 

J- X X x x 

corresponding to job We find that Z = 355. 

Return to Step 4. 


From Step 4, we find that max Z occurs at node 3. 

neNl 

Further branching takes place at node 3. The steps 4 to 7 
of the algorithm are repeated till the branching procedure is 
fathomed. A portion of the tree-diagram obtained using the 
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A portion of the Tree-Diagram obtained for Numerical 
Example 1 using Optimizing Algorithm. 
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branching scheme for this example is depicted in Fig. 2.1. 
The optimal obtained is given in Table 2.2. 


Table 2.2s Optimal Solution for Numerical 
Example 1. 


Index of Accepted 

Quantity to be 

produced 

Group 

Part 



1 

1 

30 


2 

10 


O 

2 

50 


C. 

3 

30 



1 

20 


3 

2 

60 


4 

2 

3 

7 

40 


Total 

number of units 

produced '= 

247 

Total 

time utilized 

= 

2390 mins. 


2. 4. 1.2 Model 2s Machining Speeds are Decision Variables ; 

Before structuring a mathematical model of the problem, 
we develop mathematical expressions for job production time and 
cost of production. 

Job Production Time s 

The unit production time p. . (v. .) for job J . . is a 
function of machining speed v. . for that job and is given by, 

l j 



18 


p iJ (v id 5 ’ a « + *13 + b 13 5^ 


or, 


p. .(v. .) = a. . 

^ V ij 


(i/n . - 1 ) 

■ 3 . . . b . . v. . X,J 
^id 10 11 

T7H77 

c id 


id 


¥ J. . e U 

X tj 


(2.12) 


where a. . is the preparation time, b. . is the tool replace- 

X^J X^J 

ment time, t. . is actual machining time, T. . is the tool 
X O J 

life, f\ ^ . is a machining constant, n. ^ and CL ^ are parameters 


id 


id 


for the Taylor tool life equation for job J. .. 

Let the unit production time given by (2. 12) be 

(t) 


minimum corresponding to a point v- . 

XJ 


This machining 

speed referred to as the 1 Maximum-Production-Rate Speed’ or 
! Minimum-production-Time Speedy is obtained by setting first 
derivative of p. .(v. .) with respect to v . . equal to zero. 
This yields. 


v. . 
id 


(t) 


c 13 /{b l3CV n i3- 1)} 


n. . 
id 


¥ J. . a U 
id 


(2.13: 


If w. . units of J . . are produced at a speed v. 

X^| x X 


the job production time, P. is given by, 

X J 

p i3 " S 13 +w ij Py' v 13> 


( 2 « 14 


Production Costi 


The per unit production cost, q. - (v. . ) , is also a 

X (J X vJ 

function of machining speed v. . and is expressed ass 

Xu 
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c^tvy) = a a. . + (a + (!„) + (ab.^ + Y^) ^ 


or. 


PI * ■ 

q, .(v, .) » oca, . + (a+B. ,) 


l i0 10 


10 


10 v. 


10 


+ (ab ij + *«> A ^ Vl j /n- 

id J 


(l/n -1) 


¥ J . e U 

J> J 


(2.15) 


where a is the direct labour cost, (3^ is the machining 

overhead and y . . is tool cost , respectively , 

10 

( c ) 

Let v. . x represent the machining speed correspond- 

J- J 

ing to the minimum cost point. The value of this point is 
obtained by setting the first derivative of (2.15) with 


respect to v. . equal to zero. The cutting speed, v. . 


(e) 


10 


10 


is 


called the ' Minimum-Production-Cost Machining Speed 1 and 
is expressed as follows* 

(c) _ 


V id 


( O. + p ±1 ) 

c { 3J- 


n. . 

i 13 


Cl/n ia _lKab ij +T id ) 


¥ J. . e U (2.16) 

-LJ 

Constraint s 

The various constraints for this model -are the same 
as that for Model 1. The constraint set is given by (2.1) 
and (2.5). 
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Function ; 


In this model we have two different objectives to 
consider. The first objective is to maximize the total number 


of units produced while the second objective is to minimize 


the total cost of production. These objectives are referred 
as the primary and secondary objectives, respectively. The 
expression for the total number of units produced is the same 
as in Model 1 and is given by (2.6). The expression for the 
cost of production Y,is given as; 


Y 


N N 

g g 

2 {a S.X. + 2 

1=1 1 1 j=l 


(aB . . 


+ w ij q ij( v i x ij} ( 2 - 


Optimization is carried out using the solution proce- 
dure which initially maximizes the primary objective function 
and then it minimizes the secondary objective function. 


.Analy sis ; 


Initially for all jobs J. . e U, we calculate the 

J 

minimum-production-time machining speed, v. from (2.13) 

and the unit production time is then calculated from (2.12). 

The following expression is used to calculate the job production 
time for each job considering that each unit in the lot of that 


job is produced at minimum unit production time; 

P. = s. . + 1. . p. . (v. ,( t )) 

ID ij ij ij ; 


17 


( 2 . 18 ) 



21 


Let Q t represent the total time needed to produce * 
all the .jobs to their lot sizes. Note that each unit of the 
lot is produced at the minimum-production-time machining 
speed. Mathematically, 

t % % (t)\ . 

GT = 2 (S. + 2 P. . kX; ) (2.19) 

i*l 1 j=l 10 

Depending upon the value of Q^, there are three possible 
situations. These are: Gf" = d, 0 ° > d and < d. These 
three cases are discussed separately. 

4- 4- 

i) Q = d : If the total time Q , needed tc produce all the 
units is equal to the available capacity d, of the machine, 
then all the jobs to their lot sizes are produced at the 
minimum-production-time machining speeds. Therefore, the 
optimal product mix and the cutting parameters will be given 
by the following relationships: 


. v 

V ij J 

e U (2.20) 

where U is the set of all the jobs available for processing. 


Z* = 2 1. 

J. .eU 

10 

X*. 

10 

T * 

= 1, V. . 

* 1J 

w.*. = 1. . 

10 10 


¥ J, , 
10 


4- JL. 

ii) Gt > d : When Q > d, obviously, all the units cannot be 
produced even when minimum- production-time machining speeds 
are used. Any one of the following situations may result. 


1. All the jobs are taken up for production. However, 

for one of the jobs all the units cannot be produced. 
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2. All the jobs are not taken up for production. The 

jobs taken up for production are either produced to their 
lot sizes or for one of the jobs, the desired lot size 
cannot be produced. 

The optimal product mix and machining speeds cun 
be obtained using the following approach. 

First of all, we determine the optimal machine 

loading and product mix based on the minimum-product ion- 

time speed for each job using the solution methodology 

discussed in Section 2.4. 1.1 for Model 1. Let S* and w?. 

(J-, e S*) respectively represent the set of jobs accepted 
d 

for production and the quantities of each job to be produced. 

Then, the optimal solution turns out to be the following; 

1 ¥ J. . e S* 

x* = { (2.21) 

^ 0 Otherwise 
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If slack time 6=0, i.e. if the machine is utilized 


to its full capacity, then all the jobs selected are produced 
at the minimum-product ion-time machining speeds. Obviously, 


vT . = v . (t) 

ID iD 


* J iD £ 3 


(2.23) 


If 6 >0, then we have positive slack time which, can be 
advantageously utilized to minimize the total cost of 
production. It is obvious that the minimum-production-time 
machining speed will be higher than the minimum-production-cost 
machining speed and the cost of production at the minimum- 
production-time machining speed will be higher than that at 
minimum-production-cost machining speed. So, it is reasonable 
to assume that the unit production cost is least at the minimum- 
production-cost machining speed and it increases with an 
increase- in the machining speed upto the minimum-production- 
time machining speed. Therefore, the slock time can be utilized 
in lowering down the cost of production by decreasing the 
machining speed from v. 3"^ to v. The objective is to 

determine v. < vT. < v. . 

iJ 10 - ^3 1J 


, such that the total 
cost of production is minimized. The problem can now be for- 


mulated as a nonlinear programme given below; 

N. 


S 


N. 

l 


Minimize Y = I (aS.X. + S (as . . 

1 1 j=l J 


i=l 


+ w: . q. . (v. .)) x* .} 
10 ^io 13 13 


( 2 . 26 ) 


subject to, 
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N 


N, 


i=l {SA ‘ + £ (S U + W ij P iJ (v xd» *10 1 d < 2 ' 27 > 


The inequality given by’ (2*27) represents the machine capacity 
constraint* The Equations (2. 26) and (2.27) contain certain 
constant terms which can be ignored for the purpose of 


optimization. 5 'The optimization problem can be restated as* 

N N. 

g 

Minimize Y’ * 2 S w* q, . (v. .) x* ' , 

2.-2. * =1 


(2.28) 


subject to, 

N N. 

£^ £ w* . p . . (v. . ) x* . < d 

i_l 10 10 10 10 “ < 

JU «— JL ^ — -L. 


where , 


£ 

i=l 


N. 


:Y - £ (a S . X? + £ a s . . x". . ) 

i l 10 10 


(2.29) 


(2.30) 


and 


N 


N. 

l 


d Q = d 


£ (S.X* + iT s, . x* ) 

i i 10 10 


i=l 


(2.31) 


For the problem represented by (2.28) and (2.29 ) , 
a Lagrangian function of the following form is obtained. 


1 (v,. « , p ) = 


N N. 

g l 


10 


i=l 0=1 


£ wt. q. . (v. .) x*. 


1J hLj 10 10 


N N. 


+ p, { £ w £ w* . v. ■ (v. .) x?. + k £ 
i=l j=l 10 "10 10 10 


* 


a 0 ! 


(2.32) 


where ^ and k represent the Lagrangian multiplier and slack 

variable, respectively. For optimal values of v* . , p" and k" 

XJ 

we set 9 ^d v ij = 


O. It gives, 
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or, 




-3 


3 v; - - 


7f 

10 


C*v* . 

<3 1.1 


= 0 


a+P . . 

■( —ri )Aii + (“ b ii + Y i ,M i , 


*. ( 1 /n io~ 2) 


V*. 

10 


I/nT. 


10 ’ 10 ri l0 C -/--in 

L i0 13 


* (l/n ij" 2) 


V id C. .^id 


) = 0 


From the above equation, we get, 


i ■Sf 

4’ 


p . . - ( 

K 10 n 4 „ 


v. 


1 1 ) Y • a ( fj 3 ’' 3 "") 

J ^ -i n 


l/n. . 
10 


i JL 


-i + ( U- - 1 ) b ±j ( -jM) 


vf, l/n, , 


10 


a 


(2.33) 


^ ex + q.‘* + 

3114 Vu ’ Cy 1 > (1/n id- 1)} 


n. 


13 ¥ J . . e S* 

10 

(2.34) 


Similarly, setting. 


^04 

k* 2 


O yields. 


N 


N ± 


= d - 2 E w( . p. . (v* ) x* . 

o i=1 j=1 10 ^10 10 10 


(2.35) 


and ^L/gk = 0 gives, 


l_t*k* = 0 (2.36) 

If we substitute v? . = v. in (2.33), we find that 

1J J (+) 

= 0. Similarly, if we substitute vt . = v. J in (2.33) 

J - J- J 

we obtain, p* = 00 . Hence, it is obvious that, 


4 
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M* > 0 for v7\ e (v v .^) (2.37) 

i i J U 

Since (i ' is nonnegative, it guarantees the minimization of total 
production cost. 


If \i' 


v. 


iO 


0, then from (2.36), k 5 ^ 0 and from (2.34), 

(°) « t q* (2.38) 


V ij 


¥ J. . e S 
10 


However, if \i* > 0, then from (2.36) we know that k" =0. 

The optimal values of (j,* and v* . e E. . can now "be determined 

•i J iO 

from (2.33) and (2.34), respectively, by using hit and trial 
approach. E. . is called the high-efficiency speed range. 

~L J 

Mathematical l}r , 


E. . = [v. . 

10 10 


(c) 


(t) 


iO 


3 


(2.39) 


4- -J- 

iii) Q < d s When Q < d, the optimal solution turns out to 
be the following j 


R* = 0, S* = U, X? = 1 ¥ i, Z* = Z l id> 

" x i*} ■ 1 and w Ij ■ ha * J i3 e s * < 2 - 40) 

In this situation the slack time is expressed as, 

6 = d - (2.41) 


We can minimize the total cost of production by 
utilizing the available slack time. Once the available slack 
time is obtained, the optimal machining speeds can be determined 
by formulating a nonlinear programme of the problem given by 
( 2 . 26 ) and (2.27). The optimal machining speeds will be given 
by (2.34). 
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mazing Algorithm I 


A step-by-step algorithm, for finding the optimal 
product mix and machine loading when the machining speeds 
are also decision variables, is given below s 

Step li For each nob J. . e U find v. and p./'"^ using 

~ — J ' id id - id 

(2.13) and (2.12) , respectively. From (2.18) and 
(2. 19), find . If <_ d, go to Step 2 % otherwise, 

go to Step 3. 

Step 2s The optimal solution is, 


0 = 

U 



-* 

X , = 

Id 

1 

¥ 

J id 

w * . = 
Id 

da 

¥ 

J id 

xt = 

1 

¥ 

i 


J..eS^ 

id 

If Q t = d, set v. c . = v. i* 7 ¥ J . s S * and stop. 

**“ J 3-0 1J 

For Q t < d, go to Step 9. 

Step 3 s Follow branch and bound procedure. Select an initial 

node? sot n = 1, N1 = {n}. For this node, Set R = U, 

S = 0, D = 0, s = d, Z = 0 and X ± = 1 ¥ i. Find a 

job Jw from the j cb s J. . e R, for which 
id 

(S» + s«)/l:p* + Ptp- is minimum. 


Set Z 
Step 4. 


[ d/ {(Sr + st‘t)/1tt + Prr}] and go tc 
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Find a node n from the nodes n e Nl, for which Z 


is maximum. Denote the job in n by Jrv. For the 

1 J 1 J 


node n, if Z = Z , go to Step 8j otherwise* go to 


Step 5. 


Find D = D + S7 Xr + sra~+ wrv prv 


i l 


10 "10 


where 


= max w. w. . is varied from zero to Irr till 
10 10 10 10 


D <_ d. If wo"v = 0, go to Step 6; otherwise, go to 
*** 10 


Step 7. 


Remove the job 0h-v from the list R. Set R = R-{j^v}. 

1 tj u 


If R s* 0, set Z = 2. Otherwise, find 


min { (S. X. + s. .)/!.. + pvX}, denote the corres- 
J 8 R ^ 1J lj 1J 

ponding job by J?-s and set Z = Z + [s/{ (SjXj+st-^)/!^ 

+ pf^}]- 


Return to Step 4. 


Step 7 s Branch the node n into two nodoss 


At first nodes n = n+1, Nl = Nl + {n} - {n}. 


For node n, set R = R - { Jor- } , S = S, D = D, s=s 


and Z = Z. If R = 0, set Z = Z. Otherwise, find 


min { (S . X. + s. .)/!.- + pV y } . Denote the corres- 

j li 10 10 10 

id _ 

ponding job by J-. ^ and set Z = Z + [s/{(S:rX:r + s^)/l^ 

1 J J. J. i J -L J 

+ pjV}]. Similarly at second node: n = n+1, 


id = Nl + {n}. For this node n, set R - R, 


S + {Jr^ 


D = D, s = d - D, Xr = 0 and Z = Z + wry. If R = 0, 


set Z =a Z. Otherwise, find 


d min { ( S . X, +s . - ) /l . . +p - ■ 


J . -£R 

1,1 


1 1 1J" 1J 


Cj. 
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Denote the corresponding job by Jtv and set 
z = z + [s/{(S-X, + SjjVljj + PJ 3 Ct) }]. 
Return to Step 4. 


is The optimal solution given by list n, is, 
„ * 


= S 


x. 


* 

'id 

* 


1 * J ±j e s 

■0 otherwise 


X, 


* 


Ni * 

0 if .Z, xd . = 0 

_ { a=l 10 

1 1 otherwise 


Go to Step 9. 

Find v^^ c ^ from (2.16) and from (2.12) , 

respectively. Calculate 


N 


N. 


b . * X , # (01 * . 

Q = {S A + (s 1;) + w. . p la ( >) Xy , 

If Q < d, set v* . = v < c) ¥ J. . e S* and stop. 

- 10 10 10 

Otherwise, go to Step 10. 

Step 10 : Select a job J. . e S*. Find initial value of p from 

•** 0 


the following expression, 
(3 . . - (- 1 

*10 1 

4 = — 


V. , (t) + V. . (c) 

1) T . ( jj 

n. . / ij v 2C . . 

10 10 


l/n. . 
) 10 


■1 + ( ~ -1) ^ 

n i3 10 20 ij 




) 1/n ia 


- a 


Go to Step 11. 

(p) 

Step H i For the value of p obtained in Step 10, find v. . v 

and p. . from (2.34) and (2.12), respectively. Calculate, 



30 


Q = 


N 

S 

I 

i=l 


{S.X. 

1 1 X 


(s ij + w « p ±7 } 


If Q £ d, go to Step 12; otherwise, stop. The set 
of machining speeds obtained currently is optimal. 
Step 12 ; Modify the value of ^ appropriately depending on 

whether Q < d or Q > d. If Q < d, reduce the value 
of p and if Q > d increase the value of p . 

Go tc Step 11. 

A systematic procedure for the selection of the value 
of p is presented in the form of a block-diagram in Fig. 2.2. 


Nun '. erical Example 2; 


For better understanding of the solution procedure, 
consider the single stage problem considered by Kitomi and Ham [3]. 
In all, there are ten jobs available for processing. These jobs 
arc divided into 4 groups. The basic data are given in Table 2.3. 


Solutions 


jHLs Initially, the values of v. and p. are found 

-L (J -*■ O 


for all jobs J- - e U. These are shown in Table 2.4. 

X J 


ft 


t 


5830.574 mins. Since Q > d, go to 


We obtain, Q 
Step 3 . 

Step 3? Use branch and bound procedure. 

Set n = 1, N1 = {n} = {l}. 

For node 1, we get, R = U, S = 0, D = 0, s = 3000.00 mins. 


0 and = Hi. 



Calculate value of \x from the expression given in 


Step 10 of the algorithm. Set P . = Pg = 0. 



Fig. 2.2s Block-diagram for determination of 
p and v , . . 
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Table 2.4 ; Maximum Production Rate Machining 



Speed 

and Unit 

Production Time. 

_Index 

Group 

of 

Part 

Machining 

Speed 

m/min. 

Unit Production Time 
min/ pc 


1 

223.6301 

6.715 

1 


2 

194.4333 

5.585 


1 

220.5509 

11.514 

2 

2 

151.0688 

12.667 


3 

151.4298 

7.164 


1 

144.3376 

8.219 

3 


2 

207.7381 

9.498 


1 

206.4631 

6.567 

4 

2 

188.0497 

11.912 


3 

112.1955 

8.338 


Step 4: 


Step 3 s 


S tep 7s 


Determine the values of the expression j + P 


id 


id 


¥ db^ £ R and identify the job Jj-r for which the value 
of the cxpression’TLs minimum. For the present problem, 
JV* is J, p , and the value of expression is 6.195 min/pc„ 
The maximum value of the upper bound, Z corresponds 
ro node 1. Since Z / Z, go to Step 5. 


For node 1, we find = 50 and D = 307.25 mins. 
Since w 12 > 0, go tc Step 7. 

Branch the node 1 into nodes 2 and 3? and 
return to Step 4. 


(t 
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At node 2s n = n+1 = 2, N1 = {l} + {2} - {l} = {2}. 

41* J 42» J 43* s 


Set R ^21 ? 0 22 5 '^23 5 '^31 y '^ 32 > 


S=0 } s = 3000.00 minsj D = 0 and Z = 0. Since R ^ 0 S 


determine {(S.+s- .)/l. . + p. . ^ Cji } for J. . e R. Identify 


10 ' 10 


10 


10 


the job for which the value of the expression is 

J 


J-r-r is J, and 
10 41* 


minimum. For the present problem , 
value of expression is 7. 342mir/pc. Therefore , Z = 408. 

At node 3s n = n+1 = 3, N1 = {2} + {3} = {2,3}. 

Set R = {0^27 ^21* ^22 s 23 9 ^"31 5 ^32* '^41* "^42* ^43^ * 
S = s = 2692.75 min, D = 307.25 min., Z = 50 


and X n = 0. Since R £ 0, determine {(S. + s.. .)/!._,+ p- • 


(c) 


iO' 


'10 


10 


for J. . s R. Indentify the job JV? for which the 
-*“U J 


expression carries minim-urn value. Find that Jj-r is 


Jfl for which the value of the expression is 
7.032 min/pc. Set, Z = 432. 


This time from the Step 4, we find that max Z occurs 

neNl 

at node 2. Further branching is done from node 2 following the 


steps 4 to 7 of the algorithm. Table 2.5 gives the optimal 
results obtained. 


The optimal solution shown in Table 2.5 is compared with 
the results of the same example obtained by Hitomi and Ham [8], 

It is observed that in their case, the total production amount 
was 35-6 pcs while in our case it is 372 pcs. 
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Table 2.5s Optimal Solution for the Numerical 
Example 2. 


Index 

Group 

of Accepted 

Part 

Optimal 
Machining 
Speed 
m/ min . 

Unit Pro- 
duction 
Time 
min/ pc 

Quantity to 

be produced 
pcs 


1 

173.9853 

6.81 

60 

1 


2 

134.3816 

5.53 

50 

2 

3 

126.7238 

7.14 

40 


1 

117.0901 

8.06 

30 

3 


2 

134.0509 

10.08 

72 


1 

181.4794 

6.42 

40 

4 


3 

99.0188 

7.95 

80 


Total number of units produced = 372 

Total cost of production = Rs. 20,998.87 
Total time utilized = 2,999.22 mins. 


2.4.2 Case 2s Profit Maximization s 

The problem of total profit maximization is considered 
for the following situations; 

1. The unit production time is pre specified. 

2. The machining speed which ultimately determines 
the unit production time is a decision variable. 
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3.4, 2.1 Model 3s Unit Production time is pre specified; 

The expression for job production time and the cons- 
traints for this problem turn out to be the same as for Model 1. 
The job production time is given by (2.12) while the various 
constraints are represented by (2.1) and (2.5). 

O bjective function s 


The function representing the total profit generated 
by the single stage production system under consideration can 
be written as; 


N N. 

g i 
I 2 
i=l j=l 


B. . w. . x. . 
13 13 13 


(2.42) 


where B. . is the profit per unit of job J. x. . is 0-1 
variable and is given by (2.3). The objective is to maximize 
the total profits, i.e.Z. 


Solution 



o 

* 


In this model also, we initially find Q, i.e., the 
time required to produce all the units of all the jobs. Q is 
obtained from (2.7). If Q < d, all the jobs can be produced 
to their lot sizes and the optimal solution will be; 


w. 


* 


13 

7 * 


hi * J i 3 e u 


= Z 


max 


where , 


max 


N N. 

g 1 

2 2 B.. 1. . 

i=l j=l 10 


(2.43) 
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However, if Q > d, the optimal Z < Z , Under this 

* max 

situation, we restore to B & B technique to find the optimal 
solution. 


At any stage of the B & B procedure, we select a node 

n from all the descendents for which the upper bound value 

of the objective function, Z, is maximum. Here Z is calculated 

as the sum of the lower bound value of the objective function Z, 

and max B. . [s/{(S. 1-X. + s. .)/ 1 . . + p. .}]. Note that at 
J. .eR i ij ij ij 

the fi^st node, the lower bound value Z is zero. However, at 
the subsequent nodes when a job Jra 1 is selected for produc- 
tion, Z is incremented by the quantity Brv wwv for the node 
where the selected job Jrv is assigned. The B & B procedure 
branches the node n into two nodes, one in which job J'x~ is 

X J 

in the set of selected jobs while in the other node Jrv is 

X J 

not in the list of selected jobs. For each of the descendents, 
we calculate the lower bound and the upper bound values of the 
objective function. The B & B procedure is terminated as soon 
as a descend ent with the maximum upper bound value is obtained 
for which the upper bound and the lower bound values of the 
objective function are equal. The list of jobs in this node 
gives the optimal product mix. 


Optimizing Algorithm ; 

Following is the step by step algorithm to find the 


optimal solution. 
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Step l i 

Step 2 : 


Step 3 ; 

Step 4 s 

Step 3 ? 

Step 6 ; 


Calculate Q from (2.7). If Q <, d s go to Step 2; 
otherwise, go to Step 3. 

The optimal solution iss 

_ * 


U 


w*. = 1 . . 

13 13 


* 


= 2 


max 


¥ J.. e S* 

13 

N N i 

= E S 2 B. . 1. . 

13 13 


i=l j=l 

Terminate the procedure. Follow branch and bound 
procedure. 

Select an initial nodes n = 1, Nl = {n} . 


For this node, set R = U, S = 0, D = 0, 


d. s Z 


0 


and X i = 1 ¥ i. Find a job from the jobs Jr^ £ R 

for which Yw ~ Xj + + Pj^)} 

is maximum. Set Z = Brr [d/{(S» X-*- + s-r^)/l-r~r + p-r-r} ] 
and go to Step 4. 

Find max Z and denote the corresponding node by n 

neN1 ~ „ _ 
and the job JV? in n by Jrr. If for the node n, Z = Z, 

go to Step 8; otherwise, go to Step 5. 

Find D = D + SoXr + spp' + wop; pop', where woo = max w 


13 


13 13 1 


13 


13 ' 


w. . is varied appropriately from zero to 1-yr till 

X o 1J 

/v 

D< d, If w^ = 0, go to Step 6j otherwise, go to 


Step 7. 


At node n, set R = R-{jr~}. If R = 0, set Z = Z; 

1 J 

otherwise, find a job from jobs J. . e R for which 

X j X J 
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2 = Z + B I3 [s/USjXj s I3 )/l 


ij ' 10 


Y« = max Y. .. Set 

13 J. .eR 13 
. 10 

+ Pj^lJ and return to Step 4. 

Step 7 » Branch the node n into two nodes. 

At the first node: Set n = n+l, N1 = N1 + {n} - {n}. 

Further, for the new node n, set R = R-{jr~},S = S, 

10 

s = s, D = D and Z = Z. If R = 0, set Z = Z; otherwise 

find a nob 0V-? for which Y •?--? = max Y. .. 

13 13 J. .eR 13 

10 

Set Z = Z + Bt-t [ s/ { ( St-Xt- + s^)/1t^ + p-*-*}]. 

- 10 i i 10 10 10 

At the second node: Set n = n+1, N1 = Ml + {n}. For the 


newly created node n, setR = R, S = S + {Jb"v}, s 

D = D, Xj 1 = 0 and Z = Z + wo"V Boy- 1 . If R = 0, set Z 

otherwise, find a job JV* for which Yr-r = max Y. . , 

13 13 J. .eR 13 

Set Z = Z , B-U/USjXj + Sl3 )/l I3 + P I3 }]"° 

Return to Step 4. 

St ep 8 s The optimal solution is given by node n. Set S~= S, 
Z = Z and stop. 

Numerical Example 3: 


d-D, 


Consider a single stage production system in which ten 
jobs are available for productions. The jobs are divided -into 
four groups. The basic data for the example are given in Table 2.6. 
Solution : The solution is given using the various steps of the 
algorithm: 


Q = 3105 mins. Since Q > d, go to Step 3 



Table 2.6? Basic Data Considered for Numerical Example 
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Step 3 s Follow the branch and bound procedure. 

Set n = 1, -N1 » {n} = {l}. 

For the newly cheated node 1, set 

R = { J ll s J 12’ J 21’ J 22’ J 23’ J 31 > J 32’ J 4l s J 42 ,J 43^’ 

S = 0, 2=0, s=d, D=0 and = 0 ¥ i. Calculate 

Y. . for all J. . e R. Identify the job for which Y. . 

~LJ Xj XJ 

is maximum. For the present problem Y-^ is maximum 
and equal to 0.9160 Rs./min. Therefore, is J 32 . 

Set Z = B 32 [d/{(S 3 + s 32 )/l 32 + p 32 }] = Rs. 2196. 00 
and go to Step 4. 

St ep 4: max Z occurs at node 1. Since, for node 1, Z ^ Z, 
neNl 

go to Step 3- 

S tep 3 } w 32 = 60 pcs and D = 393 mins. Since w 32 > 0, 

go to Step 7. 

Step 7 i Branch the node 1 into two nodes, 2 and 3 and return 
to Step 4. Various parameters at the nodes, 2 and 3, 
are given below s 

At Node 2s Set n = n+1 =2, N1 = {l} - {l} + {2} = {2}. 

For node n, set R = {J 11? J 21 ? J 22* J 23 > J 31 > ‘ J U1 , 

J 42 f J 43^ ! S = S, Z = 0$ s = 2400 mins, and D = 0. 

Since R £ 0, find Y . . for all J. - e R. Note that 

1 J J- J 

max Y. . is 'equal to 0.9109 Rs/min and it comes- ‘ 

J^jeR 13 

ponds to 3 ~22‘ 0 et Z = Z 3 22 £ s/ ^ S 2 X 2 ^. s 22 *) / 1 22 ^ ^22^ 
= Rs. 2178.00 



At node Set n = n+1 = 3, N1 {2} + {3} = {2,3}. 

For the newly created node n, set 

R = {J X1 » J 12 » J 2l> J 22’ J 23 5 ’ J 31» J 4l* J 42 j J 43^ * 

S = = d ~ D = 2007 mins, D = 393 mins. 

and Z = Rs. 360.00. Since R £ 0, find Y . for all 

1 J 

J-i-i e R. Identify the job for which Y. . is maximum. 
x 0 j-0 

max Y. . corresponds to job J,-, . Therefore, 

■V R 0 31 

z = Z + B 31 [s/{(S 3 X 3 + s 31 )/l 31 + P 31 }] = Rs. 2316.0 

Now we return to Step 4 and find that max Z occurs 

neNl 

at node 3. At this node, we observe that the upper 
and the lower bound values of the objective function 
are not equal. Therefore, we go to Step 5. Steps 4 to 
of the algorithm are repeated till the optimality 
criteria is satisfied. Table 2.7 gives the optimal 
solution for this example. 

Table 2.7s Optimal Solution for the Numerical 
Example 3. 


I ndex of Acc e pted 
Group Part 

1 


Quantity to be 
Produced 

Jl&cs) 

20 


Profit 

(Rs) 


340 


2 2 50 


450 


3 30 


210 


1 

20 

240 

2 

60 

360 

1 

1 

16 

3 

40 

600 


Total profit 
Total time utilized 


Rs. 2,216.00 
2,394- mins. 
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2. 4. 2. 2 Model 4i Machining speeds are decision variables;; 

For this model, the mathematical expressions for job 
production time, unit production time and production cost, 
and the constraints turn out to be the same as for Model 2. 

Equation (2. 14) represents the job production time while 
the constraints, are given by (2.1) and (2.5). The machining speeds 
corresponding to minimum-production-time and the minimum- 
production-cost are given by (2.13) and (2.16), 

Obje c tive Function , 

Our primary objective is to determine product mix and 
machining speeds such that the total profit is maximized. 
However, once the values of optimal decision variables which 
yield maximum profit are obtained, one can further manipulate 
the machining speeds such that the cost of production is minimiz 
and tho maximum profit is still realized. 

The primary objective function representing the total 
profit turns out to be the same as given by (2.42). The second- 
ary criteria of minimizing the total cost of production involves 
the consideration of total production cost expression given 
by (2.17). 

Analysis;; 

negam^m jopam ii gsa 

In this model, the analysis is almost same as done for 
Model 2. The relationships given in (2.18) through (2.41) are 
also applicable for this model. 
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Initially, Q^, is found, from (2.9) and then the time 

needed to produce all the jobs to their lot sizes at their 

maximum-production-rate machining speeds »is compared with d , 

the time available for production. Depending upon the values 
”fc 

of Q and d, one of the following situations may result . 

i) g^.dJ 

"fc 

If Q = d, then the optimal solution will be as given 

by (2.20) and the value of primary objective function Z* will 

be equal to Z 

n max 

ii ) Q^>_d.: 

If > d, then set p. . = p. for all jobs J. . £ U. 

The optimal product mix and machine loading is obtained using 
the solution procedure of Model 3. If the slack time <$ > 0, 
then this available slack time can be advantageously utilized to 
minimize the total cost of production. In view of this, the 
modified optimal machining speeds are determined using the 
solution procedure developed for Model 2. 

iii) < d : 

"j~ 

If Q < d, then the optimal machining speeds are 

determined corresponding to the minimum-production-cost. Thus, 

v*. = v. s E... The values of v* . are determined 

ID ID ID 12 

utilizing the slack time using the solution procedure des- 


cribed for Model 2. 





Following is the stepwise description of the algorithm: 

Find v. ^ and p. for all the jobs J. . e U from 
(2.13) and (2.12), respectively. Calculate Q° , from 
(2. 18) and (2.19). If condition, <_ d, is satisfied, 
go to Step 2; otherwise, go to Step 3. 

Ste]oj2 8 The optimal solution is: 


Step 3: 

mumurjgt, wUm* 



s* = u 


H - 

1 

¥ i 

if 

no = 

1 

¥ J. . e S* 

X J 

w = 
13 

1. 

1 

4 ¥ J, . s S* 

J 1 J 

Z* = 

N g N ± 

Z = I Z B. . 1. . 

maX i=l j=l ^ 


If Q t = d, set V* = V, , (t) ¥ J s S* and stop; 

otherwise, go to Step 9. 

Follow branch and bound procedure. 

Form the initial node: set n = 1, N1 = {n}. 

For the initial node n = 1, set R = U, S = 0, D = 0, 

s = d, Z = 0 and X . = 1 ¥ i. Find a job J?t from the 

1 13 

jobs J. . e R, for which Y T1 = B«/{(S T X T + s T ^)/l T? +p 

**-u J-J 1 J 1 JL lj 1 J 

is maximum. Set Z = B*--* [d/ ( (S-*- +S^)/lj-r + Py-f^^}] 

13 1 13 -13 13 

and go to Step 4. 


13 


Find a node n from the nodes n e N1 for which Z is 

maximum. Denote the job Jrr in n by J~~. If for node n, 

1 j 1 j s 

Z = Z , go to Step 8 , else go to Step 5. 
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SteJB-A* 


S tep _ 6 ; 


Step 7 i 


Step 8 ; 


Find D = D + So - * X~ + + wwy -Q/Vrv/ (t) 

i i ij ij *ij 


where 


wr~ = max w. . ? w. . is varied from zero to la-'v till 
- L J -J-J -*~J 

rv 

D <_ d. If wre* = 0, go to Step 6; otherwise, go to 
Step 7. 

In node n, set R = R - If R = 0, set Z = 2; 

10 “ 

otherwise, find a nob JV-.- for which Y = max Y. . , 

J il 13 JijeR 13 

Set Z = Z i- B b [s/^SjXj + s I3 )/l I3 + Pr^ b }]. 


-3 


Return to Step 4= 

Branch the node n into two nodes as given below,* 

At first nodes Set n = n+1, N1 = N1 + {n} - {n}. 

For the newly created node n, set R = R-{jhr}, S = S, 

i J 

s = s, D = D and Z = Z. If R = 0, set Z = Z; otherwise. 


find a job 34~? for which Y = max Y. . . 

10 J^sR J 

Set Z = Z + [s/iCSjXj + s I3 )/l^ + 

At second nodes Set n = n+1, N1 = -Nl {n} . 

For this new node n, set R = R, S = S + {Jo'c**}, 

_L J 


X^ = 0, s = d-D , D 


D and Z = Z + wr~ B~~. If R = 0, 


set Z = Z; otherwise, find a job 34*^ for which 


13 


Y i 

+ p 


max 


^ J. .eR ^ 

(4° 


Y. . . Set Z = Z + B t ^ [s/{(S t X t + s T1 )/l T - 


13 


i i 


13' 


13 


13 


}]. Return to Step 4. 
The optimal solution is, 

‘ S* = S 


x* 

X ±3 


{ 


0 


¥ J. . £ S 
13 

Otherwise 
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* 


= Z 


0 


4 - 5 


Ni 

If Z x* . = 0 ¥ i 

i=i 13 


1 Otherwise 
Go to Step 9. 

§iieJ2_9s Find - v ij^ C ^ and using (2.16) and (2.12) 

respectively. Calculate Q, 


N 


N. 

i 


Q = z" { S X" + Z (s. . + w* . p. . (c) ) x. • } 
i=l 1 1 j=i ij 13 1 1 


If Q <_ d, set v* = v. for J. . e S* and stop; 

X J 1 J X J 

otherwise, go to Step 10. 

Sten 10i Select a job J. . e S*to determine initial value of p . 

1 J 


Mathematically , 

^ -]) Y . ( 

“ij 
4 


v. , (t) + v. . (c) 1/n. . 

) 13 



a 


j Using the value of u obtained from Step 10, find 

v. and p. . for all J . . e S* from (2.34) and 
il ± il il 

(2.12), respectively. Recalculate Q from the follow- 


ing mathematical expression; 
N 

Q = 


N. 

1 


: Z { S . x ’. + Z ( s . . + w* . p . . ) X* . } . 

i=l 1 1 i=l ^ " LJ ^ 13 

If Q / d, go to Step 12; otherwise, stop. The present 

set of machining speeds is optimal. 
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Step 12 : Modify the value of \i appropriately depending on 

whether Q < d or Q > d. If Q < d, reduce the value 
of p. and if Q > d, decrease the value of p. 

Go to Step 11. 


Numerical 


Consider the numerical example 2 alongwith the 
additional information given in Table 2.8. The objective is to 
determine optimal product mix and machining speeds to maximize 
the total profit. 


Table 2.8s Basic Datas for Numerical Example 4. 


Index 

of 

Profit per unit 

Group 

Part 

Rs./pc 


1 

10.00 

1 

2 

13.00 


1 

8.00 

2 

2 

2.00 


3 

18.00 


1 

6.50 

3 

2 

8.50 


1 

11.00 

4 

2 

14.00 


3 


5.50 
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S olution s 

§L^£KLi s Determine for all jobs J. . e U, the values of v. . 

(t) ° 1J 

and p^ . These values are given in Table 2.4. 

Nov; with present- set of p Q 1 = 5830.574 mins. 

+ J 
Since Q > d, go to Step 3. 

Stgp _. 2 « Follow branch and bound procedure. 

Set n = 1, N1 = {n} = {l} . 

For node 1, set R = U, S = 0, X ± = 1 ¥ i, Z = 0, 
s = 3000 mins, and D = 0. 


Determine for all jobs Jh ^ e R, the value of Y^ 
given by expression, B i ^/{(S 1 + . 

Identify the maximum value of Y. .. For the present 

-i* J 


problem, max Y. . corresponds to job J 0 ^ and this is 
J . ,eR ^ 

equal to 2.225 Rs./min. Therefore, Z = {{^>2 ¥S 22^ Z'^2' 7 

+ = 6678.00. Now go to Step 4. 


Step 4s max Z is found in node 1. Its value is Rs. 6678.00. 
’ neNl 

Since Z ^ Z, go to Step 5. 


Step 5s D =• 323.56 mins, and w^ = 40 pcs. Since w^ > 0, 
go to Step 7. 


S tep 7 s Node 1 is branched into two nodes, viz f node 2 and 3. 

At node 2: Set n = n+1 = 2, N1 = {l} + {2} - {l} = {l}. 

Set R = R- i <J 23 > = { J p> J 12 » J 2l’ J 22> J 31* J 32* 

J 41’ ^42* J 43 ^> 8=0, s = 3000 mins, D = 0 and Z = 0. 



Since R / 0, find Y f or ' all J, . e R. Find max Y. .. 

J 1(3 J. .eR 1J 

ij 

It is observed that the maximum value of Y. • is equal 
to 1.358 Rs./min. corresponding to a job J^ 2 . Therefore, 

Z = Z + B 12 [s/Us.^ + s 12 )/l 12 + p-. 2^^ } 1 = Rs * 6344»0( 
At node 3s Set n = n+1 = 3 , N1 = {2} + {3} = {2,3}. 

For node 3, set R = R = {j^, J 21 s J 22 s J 31’ J 32 ? 

^41> ^42* J43K S = S, s = d-D = 2676.94 mins. D = D 

= 323.56 mins, and Z = 40 x 18 = Rs. 720.00, X 2 = 0, 

Since R / 0, calculate Y. . for all jobs J. • e R. 

J J 

max Y. . corresponds to job J, 0 for which 
J^eR 10 

Y 12 = B i 2 /^ S l X l * s i 2 V 1 i 2 + Pi2 (t)) ^ = 1,358 Rs */ min ' 

Therefore, Z = Z + B 12 [s/{(S 1 X x + s 12 )/l 12 + P 12 ^}] 

= Rs. 6375.00. Now return to Step 4. 


Step 4 finds node 3 with maximum upper bound value. 
Further computations are done according to the steps 4 through 12 
of the algorithm till we get the optimal machining speeds. The 
optimal solution is listed in Table 2.9. 


Table 2.9: Optimal Solution of Numerical Example 4. 


Index of 

Group Part 

Quantity 
to be 
produced 

Profit 

Rs. 

Optimal 

machining 

speed 

m/min. 

Optimal "unit 
production 
time 
min/ pc 


i 


* cOC.CO 

175.3172 

' 6.787' 

1 

2 

50 

650.00 

155.3731 

5.521 

2 

3 

40 

720.00 

127.4141 

7.121 


1 

19 

123.50 

117.7705 

8.034 

3 

2 

90 

765.00 

135.5620 

10.034 


1 

40 

444.00 

182.2416 

6.416 

4 

2 

50 

700.00 

159.3905 . 

11.225 
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Total time utilized = 2,999.66 min. 

Total profit = r s „ 3 , 998.50 

Total Cost of Production = Rs. 21, 100. 00 

2.5 HEURISTIC SOLUTION PROCEDURES: 

In Sec. 2.4 we presented algorithms for single stage 
product mix and machine loading problems considering either 
the maximization of total volume of production or maximiza- 
tion of total profit. For each of these objectives two 
variations, namely, unit production time for each job is 
prespecified and the machining speeds are decision variables, 
were considered. A solution methodology for each one of these 
situations was developed. The solution methodology utilized 
the branch and bound philosophy for generating the optimal 
solution when the total production time required to produce 
all the jobs to their desired lot sizes turned out to be grea- 
ter than the available production time on the machine. In 
order to reduce the computational effort involved in the use 
of branch and bound procedure, it was felt that some heuristic 
procedure needs to be developed for identifying the jobs which 
have every high potential for inclusion in the optimal product 
mix obtained by using the algorithms given in Sec. 2.4. The 
following heuristic procedure was developed to identify the 
set of jobs having high potential for inclusion in the final 
solution. 



As soon as the first step of any of the algorithms 
given in Section 2.4 reveals that the total available time 
on the machine is less than the total time required to produce 
all the jobs to their lot sizes, contribution of each job 
towards primary objective function, e.g., maximization of total 
amount of production or the total profit is calculated. The 
job with the highest contribution is identified and selected 
for inclusion is the list of selected parts provided the time 
required to produce the job to its lot size is less than or 
equal to the available production-time. After the first lot 
is selected, the lot with the next highest contribution is 
identified and the procedure is repeated till no job can be 
selected for inclusion in this list of selected jobs. 

If the number of selected jobs in this list turns out 
to be- less than or equal to two, this list is made empty and 
the branch and bound procedure, i.e. step 3 onwards of the 
various algorithms are followed as such. However, if the list 
of selected jobs contains more than two jobs, the last two 
jobs are removed from the list. The Step 3 of the branch and 
bound procedure is suitably modified for each of the algori- 
thms given in Section 2.4 to incorporate the heuristic 
procedure into the branch and bound approach. The remaining 
steps of the branch and bound procedures, i.e., Step 4 
onwards remain the same as in the original algorithms. The 
modified Step 3 for each case is presented below. 
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2.5.1 Caso. M _1 s Maximizati on o f Total Unit s Produced ; 

2. 5. 1.1 Case 1(a); Prespecified Unit Production Time 

Heuri_stic^ iVl^orithm ; 

The Step 3 of the original algorithm given in Sec. 

2.4. 1.1 in modified as follows! 


Step. 3..1- 


3.2: 


Step,, p., 3j 


Step 


Select an initial node: n = 1, Ml = {n} . For this 

node,, set R = U, S = 0, D = 0, s = d, 2=0 and 

X. = 1 Vi , Go to Step 3.2. 

Find from nobs J. . e R, a job J-.~r for which 

{(St Xy + st-t )/1t-» + pt-?} is minimum. Find 

1 i i 10 10 10 

'D = D + St Xt + Stt + v/tt Pt-t, where w-r- max w ,w . . 

li 10 10 10 10 10 10 

varied from zero to lr r till D <_ d. Go to Step 3.3. 

-L J 


J,.e« 


set S = S + {Jr--} 

1. . and go to Step 8. If w ?? = 0 or 


If Wtt =s Itt and D = d, oco >_> — u -r j ^ 

v J " A * < J Lj 


i0 


S + Xt =0, D = D, s 


d-D 


10 

Wt-t < 1 — r , go to Step 3.4 s otherwise, set 
10 10 

R = R - (J jjh s 

and go to Step 3.2. 

If S = 0 and w-tt = 0, then set Z*= 0, S* = 0 and stop, 

1 J 

If set S contains more than two jobs* remove the 

last two jobs selected from the list S and include 

these jobs in R, otherwise, make S empty. Set 

X. = 1 for those groups G 'of which no job is 
i ^ 

included in S. Set Z = E l- i sD={ ^ ^ s ii 

J. .eS J, xS 


N, 


+ w - p. •) + E g S.X.} and s = d-D. Find a job 
i0 10 i=i 1 1 


10 


10 


is 
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J I3 from ^ obs e R *or which {(S^ + s ij )/l ij +P iJ } 

is minimum. Set Z = 2 + [ s /{(S-rX T + s «)/l*r- + p^}3 

- 1 l i io iO 10 

and go to Step 4. 

2. 5. 1.2 Case 1(b) i Machining Speeds are Decision Variables; 
He uris tic A lgorithm : 

Step 3 remains the same as given in Sec. 2. 5 .1.1 
except the fact that instead of p. p. given by (2.18) 

is used. 

2.5.2 Ca se 2; Maximization of Total Profit ; 

2. 5. 2.1 Case 2(a) : Prespecified Unit Production Time; 

H e urist ic Al gori thm ; 

The Step 3 of the optimizing algorithm given in 
Sec. 2. 4. 2.1 is modified as follows; 

St ep 3 .1; 

Step 3.2; 


Step 3.3 ; 


Select initial list, set n = 1, N1 = {n}. 

For node 1, set R = U, S = 0, D = 0, s = d, Z = 0 


and X i = 1 ¥ i. Go to Step 3.2. 


Find a job Jr T from e R for which B-. r »/{(SjXT+s^)/lp^ 
^ u ■ 


10 


13 ' 


i i 10' 


'10 


D + S T X T + srrs 
ii i0 

is varied from 


+ p-r-*} is maximum. Calculate D 
10 

+ W I3 PJ3> where Wjj = max w ± .., w ± ^ 
zero to lrr till D < d. Go to Step 3.3. 

X J 

If = lr-r and D = d, set S = S + {04-^} 
lO 10 - L o 

2 = E B. . 1. . and go to Step 8. If w 7 rs = O or 

~ J eS « ^ 13 

w— . < 1--, go to Step 3.4; otherwise , set R = 
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3 3 + Xj - 0, D = D, s = d~D and go to 

Step 3.2. 


Step_3^4,. If S = 0 and w^ = 0, set Z*= 0, s*= 0 and stop. If 

the set S contains more than two job, remove the last 
two jobs put in list S and include them in RJ 
otherwise, set the list S = 0. X ± = 1 for those 
group G j . of which no job is included in S. 


Set Z = Z w. . B. D = { Z 

J..eS 1J J. .eS 

10 ij 


(s. . + w. . p. .) 
10 id *10 


14 


i=l 


. X. } and s = d-D. Find a jot) JV-? from the jobs 

l i u 10 


J. 


10 


e R for which B^CSj-Xj. + SjjVlj-j r Pjg} is 


maximum. 


Set Z = Z + Bj, [s/{(S T X T + s T ^)/l r , - p^}] 


l l 


10 10 


10 - 


and go to Step 4. 


2. 5. 2.1 Case 2(b ) : Machining speeds are decision variables. 
Heur istic Algor i thm s 


Step 3 remains the same as that given in Sec. 2.5.2. 1 
except the fact that instead of p. p. given by (2.18) 

1 J 1 J 

is used. 


2.5.3 Comou t at ional Experience : 

To check suitability and computational efficiency of 
the heuristic procedure, numerical examples given in Sec. 2.4. 
were solved using heuristic algorithms. It was observed that 
in all the cases, the solutions obtained by heuristic proce- 
dures and the optimizing solution procedures were the same. 



55 


Table 2.10 gives the CPU time taken on DEC-1090 computer system 
for the various examples when solved using the optimizing 
solution procedures and the heuristic approaches. It is 
observed that for all the examples considered the CPU time 
requirement was lower for the heuristic approach as compared 
to the optimizing procedure used for solving the example. 

Table 2.10; CPU Time Taken by DEC-1090 Computer 
System. 

Problem CPU^time^take^by^3^^DZlIZIIZ" 

Optimizing Solution Heuristic 

Procedure Procedure 

Numerical 

Example 1 0.13 0.10 

Numerical 

Example 2 0.16 0.14 

Numerical 

Example 3 0.13 0.11 

Numerical 

Example 4 \ 0.20 0.16 

2.5*4 An Illustration Demonstrating Use of Heuristic 
Procedure : 

Numerical example 1 presented in Sec. 2.4.1. 1 is 
considered to demonstrate the heuristic solution procedure. 
Solu tion ; 

The following is the step-wise solution of the example 
using the heuristic approach. 
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J 31* J 32 ! 


St ep 1: Q = 3105 mins, and d = 2400 mins (given). 

Since Q > d, go to Step 3. 

St on 3,»1 > Select an initial node: set n = 1, N1 = {n} = {l}. 

For node 1, set R = J 21 9 J 22 s J 2 : 

J 4l> J 42 5 J 43^’ S = 0, D = 0, s = 2400 mins, Z = 0 
and X . = 1 ¥ i, go to Step 3.2. 

Step 3. 2 s. Calculate min {(S.X. + s. .)/l. • + P M }* Its value 

is 6.55 min/pc corresponding to Jot J32* 
that D = 393.00 mins for w 32 - 60 pcs. Go to Step 3*3' 
Ste p 3 . 3: Since D < d and w^ 2 = 1^, set R = R-{J 32 ) • 


{ 3 " ^ 2 . 9 l ^" 1 0 5 


34-i 0 34 T; 9 J-z-i , J /n 9 3/, O * 3"^^} , 


'12 s °21 9 U 23 J u 31’ °41» 42* 43' 

S = S + {J 32 } = (J 32 }9 X 3 = 0 and D = D = 393.00 mms. 

Go to Step 3.2. 

Now we recalculate {(Sj_X^ + s^^)/l^j + P 4 j} ^ or 

J. . e R. It is observed that the minimum value of the 
3. 

expression is 6.8 min/pc corresponding to job J-^p’ 

We find that w^ 2 ~ 10 pcs and D — 461.00 mins. 

Go to Step 3.3. 


Since w-^ 2 = "^12 anc ^ 0 ^ 3.9 se t P 


{ J U, J 21 , J 22 ,J 2 39 


12 s 


J 32 }, X 1 = 0 and 


J 31 > J 41’ J 42 s J 43^’ S 
D = 461.00 mins. 

Similarly following Step 3.2 and 3-3 iteratively, we 
find that, R = {J 2 v 3^, J 42^ j S = ^ J ll’ J 12’ J 22 i> J 23 


J 


T t \ x=0 ¥i,D= 2289-00 mins, and 

3 2> J 32’ u 43 ; 9 ■> 


w 


42 


= 10. Since w^ 2 < 1 42 , g° to ste P 5.4. 
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‘ree-Diagram obtained for Numerical Example 1 using Heuristic 
.pproach. 
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S tep 3.4 ; We remove the last two jobs J ^ and J 43 from S and 

include them in R. Now set R = {J22’ J 31* J 41’ J 42 ,J 43^ 
S - {J n , J 12> J 22 , J 31 , J 32 }, X 4 = 1, Z * 180, 

D = 1379.00 mins, and s = 1021.00 mins. 


Find min {(S.X. + s. .)/l. - + p. J. Its value is 

J £R ^ 

0 

12.25 min/pc corresponding to job Set 

Z = 2 + [ s/{(S 3 X 3 + s 31 )/l 31 + P 31 > 3 = 263 and 
go to Step 4. 


Now we follow branch and bound procedure iterating 
steps 4 through 7 as discussed in Sec^ 2. 4. 1.1. The final 
product mix obtained is as shown in Table 2.2. Forthe purpose of 
illustration a portion the tree-diagram obtained by the use of 
heuristic algorithm is depicted in Fig*- 2. 3* 



CHAPTER III 


MULTISTAGE PROBLEM 

3.1 STATEMENT OF THE PROBLEMS 

Consider a multi-stage production system with Ns 

production stages. Each production stage has limited 

available processing time. N .types of jobs are available 

for production. Based on the part-family concept of group 

technology, the jobs have been classified into N groups. 

S 

The i-th groupj, G^ (i = 1,2,..., N ) has N. jobs. Each joh 
J . ^ (j = 1,2,..., N. ; i = 1,2,..., N_) has o lot size 

1 J 1 g 

associated with it. 

The objective is to optimally select the jobs and 
determine the quantities of each to be produced such that 
the desired criteria of optimization is satisfied subject to 
the restriction that the total processing time available on 
each stage is not exceeded. The optimization criteria is 
either to maximize the number of units produced or maximize 
the total profit. 

The unit processing time required to process a 
single unit of the job at a given stage is either pre speci- 
fied or determined based on optimal selection of machining 
speed., The optimal selection of machining speed is based on 
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minimum cost of production. 

3.2 ASSUMPTIONS: 

The various assumptions ares 

i) The jobs to be processed on a multi-stage production 
system d" not have any prespecified sequence of operations. 
Thus, job operations do net have any technological ordering. 

ii) The various jobs do not necessarily require 
processing at all stages of the production system. 

iii) Jobs do not have to be necessarily produced to 
their lot sizes. 

iv) Once a job is undertaken for production on a stage, 
all the units of this job are completed first before taking 
up the next job. 

v) Each unit of the job must be processed to completion. 

vi) Each production stage has only one machine. 

vii) Only one unit of the job can be processed on the 
machine at a time. 

viii) A limited amount of processing time is available on 
each machine (each stage). The processing time is the actual 
productive time available on the machine excluding the tine 
lost due to break-downs , etc. 

3.3 NOMENCLATURE: 

Following are the notations used for the development 
of the mathematical models: 



a ijk = P re P ara 'ti° n time for job J. . at k-th stage 
(min/pc) 

b ijk = toGl replacement time for job at k-th 
stage (min/pc). 

B i = profit per unit of job J. . (Rs./pc). 

^ijk = l~ m i a ute tool life machining speed for job 
J"ij at k-th stage (m/min) 

CUk = total utilized processing time of the machine 
at k-th stage (min) 

b k = avaa lable processing time at k-th stage (min) 

Eijk = high-efficiency machining speed range for job 
J . . at k-th stage (m/min) 

-L J 

G^ = i-th group or part-family 

i = group index (i = 1,2,..., N ) 

6 

j = job index (j = 1,2,..., N^) 




N 

N g 

h 

Ns 


P 


ijk 


= j-th job in i-th group 

= index of the production stage (k = 1,2,..., Ns) 
= lot size of job J. . (pcs) 

-i- J 

= slope constant of the Taylor tool-life equation 
for k-th production stage 
= total number of jobs available for production 
s total number of part-families or groups 
= total number of jobs in group G^ 

= total number of stages in the production system 
= unit processing time for Jh.. at k-th stage 
(min/pc) 
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(c) . . 

P±jk = rai^mum-production-cost unit processing 

(t) 


time for J. . at k-th stage (min/ pc) 

X J 


p iJ 


maximum-production-rate unit processing 
time for J. . at k-th stage (min/pc) 

1 J 


Pij(Vf^)= machining speed dependent unit processing 


S i jk 

S.. 

lk 

t ijk 

T ilk 

U 


time for J. . at k-th stage (min/pc) 

X J 

= 30b set-up time at k-th stage for J. - (min) 

X o 

= group set-up time at k-th stage for group 
G^ (min) 

= actual machining time at k-th stage for 
(min/pc) 

= tool life at k-th stage for J . . (min/ edge ) 

X 

= set of all parts available for production 


v ijk 

v i3k 


(c) 


V i3k 


(e) 


V ± 3 k 


(t) 


w. 


11 


X 


13 


X, 


Y 


{ J i/ 2 ~ N ±»* 1 - N g } 

machining speed for J.. at k-th stage (m/min) 

13 

minimum-production- cost machining speed for 
J. . at k-th stage (m/min) 

13 

high efficiency machining speed for at 

k-th stage (m/min) 

maximum-production-rate machining speed 
used at k-th stage for J. . (m/min) 

X J 

quantity of J. . to be produced (pcs) 

X tj 

0-1 type variables for 
0-1 type variable for G^ 
total cost of production (Rs. ) 
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S = set of selected jobs at any node 

= shows a candidate for acceptance or rejection 
[A] = is a Gaussian notation implying a greatest 
integer less than or equal to A. 


3.4 MATHEMATICAL MODELS AND SOLUTION PROCEDURES : 

Mathematical models for the multistage problem stated 
in Sec. 3*1 are presented along with optimal solution methodo- 
logies. The maximization of profit is considered as the cri- 
teria for optimization. Separate models and solution methodo- 
. logics are presented for the following two cases : 

i) The unit processing time is prespecified. 

ii) The machining speed which ultimately determines 
the unit processing time is a decision variable. 

3.4.1 Total Profit Maximization When Unit Processing Time is 
Pre specified : 

3. 4. 1.1 Constraints i 

The following constraints are considered: 


i) Machl ne utilization constraint ; 

This constraint ensures that the total time needed 
to process the various jobs at any stage does not exceed the 
available capacity of that production stage. Mathematically, It 

can be written as: 


g 

£ 

i=l 


(S., X. + _ 
lk l 0 


N. 

£ 


JL p 


. x. . ) < 

x ijk “ 

k = 1, 2, • • 


Ns 


(3.1) 
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whore , 


^ijk •3°^ Processing time for processing w^.. 


units of job J . . at the k-th stage and is given by, 

X J 




(3.2) 


X ij 


jk ~ G ijk + w i j p ijk J 

1 if job J^. is selected for production 


{ 


and. 


13 

0 if job is not produced 


N-s 


if I x. - > 1 


h 


3=1 

N,- 


13 


(3.3) 


(3.4) 


0 


if E x. . = 0 
j=l 10 


ii) Lot Size Constraint ; 


The quantity of job produced should not 


exceed the lot size specified for it. Thus, 


0 w ij ^ hi 


J" . . £ U 
13 


(3.5) 


3. 4. 1.2 Ob.iective function 


The total profit can be expressed as£ 
Z 


l Z B. - w X 

i=l j=l 13 13 3 


(3.6) 


where, B.^ is the profit per unit of job >' i3 is 3 0-1 

variable defined by (3.3). The objective is to maximize the 

total profit, i.e. ,Z. 
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3. 4.1. 3 Solution procedure ; 


Lot represents the time needed, to process all 
the jobs to their lot sizes at the k-th stage of production. 
Then, 


'4 


k 


N g N ± 

is]. ^ SjLk + j=i ' Sj -3k + ^j p ijk^ 
k = 1, 2, . . . , Ns 


(3.7) 


If for all the production stages the condition 

jr ]/• 

Q d is satisfied, then all the jobs will be produced to 
their lot sizes. Thus, the optimal solution will be: 


where, 


* 

z 


= i 

= l 

= z. 


Ij 


¥ J. . e U 
id 

¥ J. . e U 
10 


max 


N 


(3.8) 

(3.9) 

( 3 . 10 ) 


max 


i=l j=l 


N. 

l 


If the condition Q k < d k is not satisfied for any 
of the production stages, then any one of the following 
situations may arises 

i) All the jobs except one of them, are produced to 


their lot sizes. 

ii) All the jobs are not selected for production. 

In this situation two cases are possible: 
a) Jobs selected for production are produced to 


their lot size. 
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t>) One of the selected jobs can not be produced 
to its lot size. 

In all the above situations, it is obvious that 
^ ^ ^max* ^ branch bound (B & B) procedure is developed 

to determine its optimal product mix and machine loading for 
each of these cases. 


The B & B procedure is used to branch a node n for 
which the upper bound value Z, is maximum. For each of the 
descendant obtained after branching a node, the upper bound 
value Z and the lower bound value Z is determined. To obtain 
Z, a job from the remaining jobs is identified for which 


B. ■ N. . is maximum. The corresponding job is denoted as JVv. 

X J X J 1 J 

Mathematically, N. . is expressed ass 



min [ 
k 



+ p iok^ 


( 3 . 11 ) 


where s^ is the time available for production at that node. 
The job Jro' is included in the list of selected jobs in one 

X J 

of the two nodes obtained from branching of a node having 


maximum upper bound value. For the descendent, for which the 

job Jrr is included in the list of selected jobs, the lower 
i<3 

bound value Z, which at the beginning of B & B procedure was 
set equal to zero, is incremented by Bro- wo-r. The list of the 
selected parts corresponding to the node for which the upper 
bound value is maximum among all the nodes and equal to the 
lower bound value at the same node, gives the optimal product 
mix. At this point the B & B procedure is fathomed. 



3 . 4. 1.4 Optimi zing algorithm : 

The various steps of the algorithm ares 
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Step 1 

Step 2 

Step 3 


Step 4 

Step 5 


Find Q k , for each stage from (3.7). If \ < d k ¥ k, 
go to Step 2; otherwise, go to Step 3 , 

Produce all the jobs to their 1 lot sizes and terminate 
the procedure. The optimal solution is given by the 
f :1 lowing relationships s 


x 


i 3 


1 * 
1, W. . 

10 


1 . . ¥ J e U and Z ' » Z 

lj ij max 


Follow branch and bound procedures 

Select an initial node? set n = 1 , N 1 = {n}. 

For this node, sot R = U, S = 0, s k = d k 

Find, 
s. 


D k a 0 ¥ k, Z 


0 and X. = 1 ¥ i, 


Find N . = min [ r —f"; _ 

10 k fik X i + s i.ik 

i±r 


] ¥ J. e R- 


+ p 


ijk 


Bind max g p , . and denote the corresponding job 
J ij eR 13 10 

by Jrr. Set, Z = Bt” Ntt and go to Step 4. 
i j 1313 

Find a node from the nodes n e N1 for which Z is 
maximum. Denote this node by n and in n by 

For the node n, if Z = Z go to Step 8 ? otherwise, go 
to Step 5. 

Find for each stage D k - D k + Sj; k Xy + s£j k + W 33 p^. 


max w. w. . is varied from zero to Ir'r 
1 J 1 J 


where woo- 
J 

till D k < d k ¥ k. If w^v = 0, go to Step 6 ; otherwise, 
go to Step 7. 
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Ste p 6 ; 


St ep 7 s 


Step 8 s 


Remove the job Jrv from the set R of the node n. 

Set R = R - If R = 0, set 2 = Z; otherwise, 

find, j ma ^ ®fj ^j_j* denote the corresponding job 

ky and set 2 = 2 + Bj-j Njj. Return to Step 4. 

Branch the node n into two nodes, as given belcws 
At first Nodes 


Set n = n+1, N1 — N1 + {n} — {n}. 

For the newly ere ted node n, set R = R - { , 


S = 0, Sy. - s k and D k = D k 
and Z = Z. If R = 0, set Z 


¥ k, N. . = II. . ¥ J. . e R 
U ij 10 

= Z* otherwise, find 


max B. . N. denote the corresponding job by J-t 
. eR lj *3 

and set Z •- Z + Btt IU. 

10 10 

At second Nodes 

Set n = n+1 and N1 = N1 + {n}. 

For the newly created node n, set R = R, S = St 

±. J 


D k = D k and s k “ d k ~ D k ¥ kj 


X. =0 and Z = 2 + 


+ B'ro' wre', If R = 0, set Z = Z$ otherwise, recalcu- 
10 10 

late the value of N. . ¥ J . . e R and find max B N . . 

xj ij ij ij 

Denote the corresponding job by and 

S6“t Z = Z + . 

Return to Step 4. 

The optimal solution is' contained in node n. 


* 


Set, S = S, Z 


Z and stop . 


The value of wo"v required in Step 5 is determined using 
X 3 

the following relationship. 
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= 

min [ 
k 

ikZJk. 

- Sr, 

XK 

Ij 


P ijk 


If, wrr 
U 

\ 1 

ij 9 

we set 

il 


3. 4. 1.5 Numerical Example 5 : 

Consider a two stage production system for which ten 
types of jobs arc available for production. These jobs are classi- 
fied into lour groups. The basic data for the problem is given 
in Table 3*1. The objective is to determine optimal machine 
loading and product mix. 

3 « 4 . 1 . 6 Soluti on s 


Using the various steps of the algorithm, we obtain, 

Stop 1 : = 3105.00 mins, and = 2995.00 mins. 

Since ¥ k, go to Step 3. 

St ep 2 s Select an initial node? Set n = 1, N1 = {n} = {l}. 

I 1 on node 1 , set R = ^2.2* ^2.1.* ^22 * ^23* ^31 s ^32’ 

^ 41 » ^42’ J 43 ^ » S = 0, s 1 = 2400.00 mins, s^ = 2400.00 mins, 
Z = 0 and X. = 1 ¥ i. Calculate N. . ¥ J. . e R. Find 

X X J IQ 

max B. . N. .. Its value is Rs. 4080.00 corresponding to 

J. .ell 1J 

1 J 

job J-^ 0 . Therefore, Z = B-^ N-^ = Rs. 4080.00. 

Go to Step 4. 

Step 4 ; Since for node 1, Z / Z, go to Step 5. 

Step 5: For w-^ = 10 pcs, = 68.00 mins and = 94.00 mins. 

Since < d k ¥ k, go to Step 7. 

Step 7 : Branch node 1 into two nodes, viz, nodes 2 and 3. 

At node 2: Set n = n+1 =2, N1 = {l} + {2} - {l} = {2}. 
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Tabic 3.2s Optimal Solution of Numerical Example 5. 


Index of 

Group 

Accepted 

Job 

Quantity to 
’be produced 

Profit 

obtained 

Rs . 


1 

30 

195.00 

1 





2 

10 

160.00 


1 

20 

390.00 

2 

2 

50 

625 . 00 


3 

30 

270.00 

3 

2 

60 

840.00 

4 

3 

40 

592.00 

Total 

number of 

units produced = 240 

Total profit obtained 

= Rs. 3,072.00 

Total 

processing time utilized: 


At stage 1 

= 2,394.00 mins. 


At stage 2 

= 2,386.00 mins. 
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and '6-0, We find that .s^ = 82 = 2400,00 mins, 

Z = 0 and = 0 ¥ i. Since R £ 0, find 

max B. . N. .. Its value is Rs. 3025.00 correspond- 
J iO SR 10 ° 

ing to a job J 22> Set Z = Rs. 3025.00. 

At node 3s n = n+1 = 3, N1 = {2} + {3} = {2,3}. 


R 


S 


^^11 J ° T 21 > ^ 22 ’ ^23’ ^”31’ ^ 32 J ^4l» ^42» ^*43^ 
{J 12 }. We find that s^ = 2332.00 mins, 


and 


s 2 = 2306,00 min, D-^ = 68.00 mins,, D 2 = 94.00 mins, 

Z = Rs. 160.00 and X x = 0. Recalculate the value of 

i ^ J . i £ R. Find max -g _ -p its va lue eorres- 
XJ Xd J id eR 

ponds to a jo* 3 J 22 . ihenei 0336 ? ^ = Z + B 22 k 22 = 


- Rs. 3110.0. 


Return to Step 4. 

From Step 4, we find that max Z occurs at node 3. 

nshl 

Therefore, the further branching takes place at this node. 

The steps 4 to 7 of the algorithm are repeated till the 
branching procedure is fathomed. The optimal solution obtained 

is given in Table 3.2. 

3.4.2 Total Profit Maximization wh en Machining Speeds 
are Decision Variables ? 

Before developing a mathematical model and solution 
procedure for the problem, we develop mathematical expressions 
for O'* processing time and cost for processing a single unit 


at a given stage of production. 



3.4.2,] Job processing tim e; 

The unit processing time, p^ k required to 
process a single unit of the 30 b at the k-th production 
stage is a function of machining speed v ± ^ k for that 3 ob at 
k-th production stage. Mathematically « 


•'ijk ^ V iQk^ 


, . ji.ik 

a. + t. ., + b. . v sfr- 1 *— 

13k 13k 13k T i3k 


or. 


/ l/n . -1) 

' n u 


V J-/ li - ‘T, 

.lull — bl 

a 3 K Ciik « k 

(5.12) 


“ijk 

V J e V 


whom, a iSk , b iJk and t.. k , respectively represent the 
propurnti h time, tool replacement time, and actual machining 
time on li-th production stage. T^, r- i<jk and C iJk are the 
parameters of the Taylor tool life equation. 

Let v. . represent the machining speed at which 

p . .. (v., k > is" 1 ” minimum, v^^ is obtained by setting the 

first derivative of P ljk (v 1Jk > with respect to v. Jk equal 
to zero, 


(t) _ r /U ) n± j i = 1 . 2 , 

v ijk “ ¥ n ijlc 


mono iV 


(3.13) 


» % 


3 ~ 1 ? 2 $ * « * 
k = 


~ ■ -u j are processed on k-th stage at a speed 
If w . . units of 30 b Jy are 

2-0 ** 


v 


Ijk 


the job processing time, Ip k . i® given by, 


X > 

13k 


s ijk + W ij P « k (Vi 3 k> 


(3.14) 
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*3. 4. 2. 2 Production Cost: 


The cost of processing one unit of job J. . at stage k, 
^i jk v i jk^* a ^- s0 a function of the machining speed v. ... 

X J K 

at that stage. Mathematically, it is expressed as; 

q ijk (v iJk > = % a ldk + (o k + P ±ak ) 

, ^ 11k v 11k (1 /n ldk -l) 

+ ( “k b ijk + U d k) -f-fc 

ijk 13K 

(3.15) 

where (3^..^, Yijk ? respectively represent the direct labour 

cost, machining overhead cost and tool cost for job J. . at 

i<D 

production stage k. 

( c ) 

Let v^^ represent the 'Minimum- Production-Cost 
Machining Speed 1 . It is obtained by setting the first deri- 
vative of (3.15) with respect to v^^ equal to zero. Mathe- 
matically, 


v ijk 


C,.,„ ( 


(a k + h.ik ) j 1 

_ ly W hjk + Ujk^ 

(3.16) 


3. 4, 2. 3 Objectr 


Our primary objective is to determine product mix 
and machining speeds such that the total profit is maximized. 
However, once the values of optimal decision variables which 
yield maximum profit are obtained, the machining speeds can 
be further manipulated such that uhe cost of production is 
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l 


minimized without sacrificing the maximum profit. The primary 
objective function turns out to be the same as given by (3.6). 
The secondary criteria of minimizing the total cast of produc- 
tion involves the total production cost expression which is 
given as: 


Y = 


Ns N g 
E 

k=l i 


N. 

1 


El {a k S ik X i - Z Ca k s lak 


10 H ijk v ijk 


^ V iik^ X ijt^ 


(3.17) 


3 . ^ . 2 . 4 Analysis 


I* or all the jobs, the minimum-producti on-time machining 
speed and unit processing time at the various production stages 
are calculated from (3.13) and (3.12) .respectively. The follow- 
ing expression is used to calculate the job processing time for 
each job at each stage considering the fact that the each unit 
in the lot associated with that job is produced at minimum unit 
processing time. 


p (b) _ g ' + 1 o s (v (b)'s 

F ijk s ijk lj p ijk tv ijk ' 


( 3 . 18 ) 


X 

Mathematically, Q k , the time needed by k-th production stage 
to process all the jobs to their lot sizes at minimum-production- 
time machining speeds, can be expressed as, 

N 


N. 

< - £ <hk * £ hok 


(t) 


), k = 1,2,..., Ns (3.19) 


Depending upon the value of Q k , there are three 
possible situations s Q k " fc = d k ¥• k, > d k for at least 
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one of the production stages and Q, t < d, ¥ k. These three 
cases are discussed separately. 

1} Q K. = .. d k.t k: 

= d^ ¥ k (k = 1, 2, » . . , Ns), then produce all 
the jobs to their lot sizes (at the minimum-producti on-time 
machining speeds). In this situation, the optimal product 
mix and the cutting parameters are given by the following 
relationships; 

X ij 1 

w. .* = 1. . f ¥ J, . e U (3.20) 

10 10 j 10 

and "^"ijk = "^"ijk ^ = 4 ■ ‘ 5 ^ 

ii) > d k for at least one of the prod uc tion stages ; 

In this case, the processing time available on at least 
one of the production stages is not sufficient to process all 
the jobs to their lot sizes. In fact, any one of the following 
situations may result: 

i) All the jobs, except one are produced to their lot sizes. 

ii) All the jobs are not produced and 

a) the jobs selected are produced to their lot sizes. 

b) all these jobs except one are produced to their 
lot sizes. 

In all the above situations, the optimal product mix 

be obtained using the following approach. 


and machining speeds can 
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First of all, we determine the optimal machine loading 
and product mix based on the minimum-production-time speeds at 
each stage using the solution methodology discussed in Sec. 3.4.1. 
Lat and (j\ a s ) respectively represent the set of 

jobs accepted for production and the quantities of each job to 
be produced. Following is the optimal solution. 


x.'. 

10 


{ 


Further, 


X i 


* 


- { 


1 ¥ J. . £ S* 

10 

0 Otherwise 


% 

1 If 2 x. . > 1 

0=1 10 ~ 

0 Otherwise 


( 3 . 21 ) 


( 3 . 22 ) 


per each production stage, the total time consumed Cu k and 

the slack time 6^ are calculated from the following expressions, 


Cu, 


N N i 

E'{S ik x; + ^ (s ljk+ w*. p.. k 


and 


i=l 

d k 


(t) 


> "iO > 


Cu. 


k 


(3.23) 

(3.24) 


If by = 0, then all the selected jobs are processed at their 
minimum-production-timo machining speeds (at that stage). 
Obviously, 

(t) it -T c- tt (3.25) 


* 


= v. 


¥ J. . 6 U 
IQ 


ijk v ijk 

If 6 k > 0, it implies that there is some slack time 
available at the k-th stage. This slack time can be advantage- 
ously utilized to minimize the total cost of production. Since 


v ijk 


(t) 


> v (c) and q.. k (v.^O > d ijk (v^'), it is 


(tx 


(c) 
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« 

reasonable to assume that the cost of processing a single 

unit of a job is least at the minimum-production-cost machining 

speed and. it increases with an increase in the machining speed 

upt 1 ' 'the minimum-production-time machining speed. Therefore, 

the slack time can be utilized in lowering down the cost of 

( is ) 

production by decreasing the machining speed from “t 0 

(c) • * „ (t) 


V, 


Thus, the objective is to determine , 


v ijk 


such that the cost of production is minimized. 


ijk 

< v.* 

1 , J K. 

The problem is formulated as a nonlinear program as given below; 

N 


< v .>> 
10 


N. 

Ns g * 1 , 

Minimize Y = I 2 {a k + _2 s j_jp 

k=l i=l 0=1 


+ w A % jk > nj > 


( 3 . 26 ) 


subject to 
N 


g 


A {5 ik 

1=1 


* 


N i 

+ ^ ^ s ijk * w io 


?ijk ^ v ijk^ ; x ij ^ 


Ns 


(3.27) 


k — 1 5 2 9 • « • . 

There are certain constant terms in (3*26) and (3*27) 
which can be ignored for the purpose of optimization. The 
problem can be restated ass 


Minimize Y 1 


» r N N. 

Ns g 1 , v * (-2 

S. A w i3 q idk (v i d k) x ijk (3 
k=l i=l 0=1 


subjcct to, 
N 


N. 


i=l d 5 * 1 


* 

w ij p 


4 


ijk ^ v i jk^ x ij 


< do,, k = 1,2, . 


.28) 


Ns 

(3.29) 
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where, 

Y’ 


Ns N g 

= Y " k=l ifi S lk h* ^ £ % s iJk x* a ) (3.30) 


and . 


N N. 

1 f * _o f * 1 

k = d, *- Z (S..X.+Z s. x * \ v -to M „ t-x 

k i=l 2 -;_i i,ik i .1 * & l*2 9 * e , ? Ns (3*31) 


d=l “iOk -ij 


For the problem represented by (3.28) and (3.29), a Lagran- 
gian function of the following form is obtained. 


L (v iok- V 


Ns Ns 

= 2 E Z w ’ a (v W* ' 

k=l i=l j=l x d q iJ k 1 ijk' x ljk 

Ns N g N. 

4* 2 U n { E E w . . p . . . ( v . . ) x J . 
k=l 1 i=i j=l i jk ij 

+ K k 2 - do } 

whore n k and K k represent the Lagrangian multiplier and slack 
variable, respectively. For optimal values of v ± * k , & K* , 

we set 32 / 0 v ijk = 0. It gives 


15j-i k ._ ( liak.2 + l_ P i,jk (v i!ik l 

- 3 V ^ 


0 


ijk 


or. 


c 


a k + P-i 


^ >/\j k + (« k b 13k *V 1)} J 


v iak 

* v ». (1/n i 3 k- 2) 

0 

r 1 J 

Id 
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* f- , iyk b i.1k v i.i 


+ LI f_ '11, . 

K k ' 777 * + 


m_4dk 
1/n. ., 
p ijk 

ijk 


* ( 1//n iQk" 2 5 


) = 0 


From the above equation, we get 


, , 1/n. .. 

(1/n. . . -1) y . ( v * /C ) 

7 ink j r i,j kV i,ik 7 i,ik' 


- 1 + b. , k (v*. k /c k ) o* 


1,2,..., Ns (3.32) 


V ijk 


C. .. { 
ijk 1 


a k + u k + ^ 


: n iik 


^ a k + ^k" 13 ijk + Y iok^ 1//n i,jk" ;L ^ 


i = 1,2, • . • , N 


3 53 1 $ 2. 9 • o • y 

N ^ 1 s 2 9 <r 9 * } Ns 


(3.33) 


Si milarly 


» setting 3Z/^]i k = 0 yields. 


iM ^ r 

g N. 

do k ' 3l •/ W « P « k <V ijk ) k*. , k = 1,2,..., No 

(3.34) 


and 0 Z/0K, = 0 gives. 


M' — 0 -y k = JL j 2j « * « j Ns 


:3.35) 


If we substitute v* . = v. ,/- c ^ in (3.32), we find that 

ijk ijk 

q’ k = 0. Similarly, if we substitute vf_^ = v^ k 7 in (3.32), 
we obtain \xf = 00 . Hence, it is obvious that, 


h k > 0 for W. k e (v ijk (c) , v 13k W), k = 1,2..., Ns 


(3.36) 
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Since is non-negative, it guarantees the minimization of 

cost of processing selected jobs on k-th production stage. 

If = 0, then from ( 3 . 33 ), v* jk = v.. k {c) ¥ j e S* . 

However, if p^ > 0, then the optimal values of and 

v i"jk e ^ijk ^ or eac ^ °^’ pnoduc.ti.on stage is determined 

from (3.32) and (3*33) > respectively, by using hit and trial 

approach. E. is the high-efficiency machining speed range for 

job J„. . at stage k. Mathematically, 

1 0 


E. = [v. .. ^ 
13k L 13k 


(t) 1 
V. M j 
i3^ 


(3.57) 


lii) 0* < ¥ k : 


When Q, t < d. t ¥ k, the optimal solution turns out to be 
K " w " K 


the foil owing j 

H ' = > , S x '= U, Xf = 1 ¥ i, ‘x% = 1 and w* 


no 


1. . ¥ J . . £ S' 
iO no 


(3.38) 


For the production stages for which ^ - d^ , the optimal 


me chining speeds vl k are 




v ijk = V i0k 


¥ J • • a S'*, 
i 0 


However, for the production stages for which < d^ * 
the optimal machining are determined as follows: 

We calculate the slack time, 


6 k = d k " Q k 
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The optimal machining speeds are determined by formulating 
a nonlinear program of the problem a,s given by (3.26) and 

(3.27). Hit and trial approach is used to obtain the optimal 
solution from (3.32) and (3.33). 

3 * 4 , 2 ♦ 5 Opt imi zing Algorith m : 

The following is the step by step description of the 
optimizing algorithm: 

Step.,!.: Find Q k t > for each stage from (3.19). If Q k t < d^ ¥ k, 
to Step 2, otherwise, go to Step 3. 

St op 2 : Tho optimal, solution is: 




Find max B. . N._ and denote the corresponding job by 

1 j 1 J 


J-r-r 

10 


J. . £ R 
iO 

Set, Z = Btt N«r~r and go to Step 4. 
10 10 


Find a node from the nodes n £ N1 for which Z is 
maximum. Denote this node b 3 ^ n and J-£^ in n by JArp 
For the node n, if Z = Z, go to Step 85 otherwise, go to 
Step 5. 

Step, pjt Find for each stage D k = D k + Xv + srv R + wer pej, 

max vn w. . is varied from zero to Ivv 


where wo*v - 

.L J 


10 s 10 


till D k < d v ¥ k. If wrv = O s go to Step 6 ? otherwise; 


A k 


10 


go to Step 7. 


Step C; Remove the job Jvo' from the set R of the node n. 

Set R = R - If R = 0, set Z = Zj otherwise, 

X0 

find max B N . „denote the corresponding job by 
OhjCR ij i 0 

Jrr and set Z = Z + Brr Nr*r. 

13 10 10 

Return to Step 4. 

Stop 7‘* Branch the node n into two nodes, as given below. 

At first nodes Set n = n-;-l, N1 - N1 + {n}”{n}» 

For the newly created node n, set R — F— S 0s 


s k 


s k and D k = \ * k > N. . 


10' 

N ¥■ J. . £ R and Z = Z. 
ij 10 

N, 


If R = 0, set Z = Z; otherwise, find ^ max 
denote the corresponding job by Jjg ^nd set 
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At .second nodes Set n = n+l s N1 = N1 + {n}. 

tor the node n, set R = R, S = s = D^. and 

s k “ P k ^ X T ~ 0 and Z = Z + B-xy N'T'y. If 

R = 09 set Z = Z; otherwise, recalculate the value of 

R-i i ¥ J. e R and find max B. . N. .. Denote the 

J 1J J. . e R ^ ^ 

, . il 

corresponding job by j— and set Z = Z Br-r N¥-». 

id - il il 

Return to Step 4. 


Stpj)_ n> 8_: 


The optimal solution is contained in node n. 
* 


3tpp m _ 9 i 


Set S = S, Z A = Z, Q. 
go to Step 9. 

If a 


¥■ k, k = 1 and 


“k 


d,., set v* dk = v.^^ V e S * and go to 


A k’ 


(c) 


Step 13; otherwise, if Q k < d, 9 find v. - k w and 

K K 1JK 

p^ k from (3.13) and (3.12) respectively. Calculate 

(c) 


Qi 


k 


N g 

Z 

i=l 


V. 

N. 

1 

s ik x i 

+ ^ (s ilk + W ij p : 

ijk 

v. ^ ¥ J. . £ S* 

ilk il 


x 


*. } 
il J 


w k -• k J 

Step 13; otherwise go to Step 10, 

Step_10, Select a job J. . e S *, Find initial value of \x , from the 

-L J 

following expressions 




(t) + v (c) l/n iik 

««, - <■/-«« -« >«, < 


V. 


-1 +(l/n^ j k ~l) b i j k ( 


2C iik 


-) 


Go to Step 11. 
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For the value of p, obtained in Step 10, find v„. 

X IV 

an 4 P-; from (3*33) and ( 3 * 12 ), respectively. 


(e) 


Calculate , 


Q, 


l k 


Ng ^ N ± 

ifl ^ S ik X i + ^ ^ s i jk + w ij p ijk^ x ij ^ 


If Q,_ = d n _„ set v*., = v. .. 


k ~ ^k s ocu v ijk = v ijk 
Step 13 5 otherwise, go to Step 12. 


e S * and go to 


Step 12 ; Modify the value of p^. appropriately depending on 
whether Q^. < d^ or > d^. If < d^, reduce the 
value of p^. and if > d^., increase the value of p^.. 
If k < Ns, set k = k+1 and go to Step 9? otherwise, 
stop. 


3 * 4 . 2 . 6 Numerical 


Consider a two stage production system for which ten 
types of jobs are available for production. These jobs are classi- 
fied into four groups. The objective is to determine optimal 
machine loading and product mix when the machining speeds are 
decision variables* The basic data for the example considered 
is given in Tables 3*3(a) and 3 . 3 (b). 


3. 4. 2.7 Solution; 

To find the optimal solution of the problem given in 
Sec. 3. 4. 2. 6 , we follow the various steps of the optimizing 
algorithm given in Sec. 3. 4. 2. 5. The optimal solution is presented 
in Table 3.4. 
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Table 3.4: Optimal Solution of Numerical Example 6. 


Index 

of Accepted 

Quantity to be 
produced pcs 

Machining 

m/min 

s^eed 

Group 

Job 


At Stage 1 

At Stage 2, 


1 


8 

143.106 

156.739 

1 







2 


50 

133.109 

115.709 


1 


100 

127.29 

179.907 

2 

3 


40 

110.596 

131.009 


1 


40 

163.117 

157.462 

4 







2 


50 

138.144 

171.259 



Total 

Profit = 

= Rs. 3,390.00 



Total 

Production Cost 

= Rs. 33,879.48 



Total 

Processing Time 

Utilized: 





At stage 1 

= 3,000.00 

mins* 




At stage 2 

= 2 , 500.00 

minso 


3.5 HEURISTIC SOLUTION PROCEDURES: 

In Sec. 3.4, solution methodologies for two variations 
of the problem, namely, prespecified unit processing time and 
machining speeds are decisions variables, were developed. The 
objective was to maximize the total profit. These solution 
methodologies used the basic concept of B & B for providing 
the optimal solution when at least on one of the production 
stages, the total processing time required to process all the 
jobs to their lot sizes was greater than the available processing 


time on that stage. 
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It was realized on solving many multistage production 
system problems using the solution procedures given in Sec*3.4, 
that the computational effort increases with the number of 
the production stages and the number of jobs available for 
production. Therefore, an effort was made to reduce the compu- 
tational burden involved in the use of B & B procedure. This 
was achieved by reducing the size of the problem to be handled 
by B & B procedure. The reduction in the size of the problem is 
achieved by identifying a set of jobs which seem to have high 
potential for appearing in the optimal solution generated by 
the solution methodology developed in Sec. 3. 4. 1.3. The heuris- 
tic procedure is described as follows. 

As soon as, the first step of the algorithms, described 
in Sec. 3.4, indicates that there is at least one production 
stage on which the available processing time is less than the 
total time required to process all the jobs to their lot sizes, 
contribution of each job towards the objective function, i.e., 
maximization of total profit is calculated. A job with the 
highest contribution is identified and is selected for inclu- 
sion in the list of selected jobs provided the available processing 
time is sufficient to process this job to its lot size. After 
the first job is selected, the next job with highest contri- 
bution is identified and the procedure is repeated till no job 
can be selected for inclusion in the list of selected jobs. 
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If the number of selected jobs in this list is less than 
three, the list is made empty and the B & B procedure, i,e.. 
Step 3 onwards of the algorithm given in Sec. 3.4, are followed 

I* 

as such. However, if the list of selected jobs contains more 
than two jobs, the last two jobs are removed from the list. 

The Step 3 of B & B procedure is suitably modified for each of 
the algorithms given in Sec. 3.4 to incorporate the heuristic 
procedure into the B & B procedure. The remaining steps of the 
algorithms, i.e.. Step 4 onwards remain the same. The following 
section gives the modified Step 3 for the various cases. 


3.5.1 T otal Profit Maximization 'When Unit Processing Time is 

Pre specified ; ~ ~~ 

3. 5. 1.1 Heuristic Algorithm i 

The Step 3 of the original algorithm given in Sec. 3. 4. 1.4 


is modified as follows: 

Step 3.3f t Select an initial node: n = 1, N1 = {n}. For this node. 


set R = U, S = 0, s^ « d k and D R 
X. = 1 ¥ i. Go to Step 3.2. 



Calculate, 


0 ¥ k, Z = 0 and 


N ij = ”f n [s k /((s ik x i + s l3k > / 1 i3 + p ijk^ ¥ J ij E R 

Find max B. . N. . and denote the corresponding job 
J. xR 10 1J 

by Jtt. For each production stage, calculate 
i J 

\ - D k + s lk x i + s l3k + W I3 p l3k’ where W I3 = nax 
w. . is varied from zero to l— till 1 < tl » k. 

JL J -*“U 

Go to Step 3.3* 
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If w T -3 = 1 7 , and D. « 


"ij = ij and D k = d k ¥ k » set s = s+ { J ij}* 

T Z * B ii hi and S° to Step 8 * If wm = 0 or 
Jf^eS XJ x 0 10 


W I3 < 1 I3 5 so to step 3.4$ otherwise, 

set R = R _ {j„ }, s = S . {Jjj}, = B k>Sk - d k - D] 

and Xj = 0. 

Go to Step 3.2. 

§Ji§IL 3,^4: If S = 0 and = OJ set Z*= 0, S* = 0 and stop. 

If set S contains mor@ than two jobs* remove the 
last two jobs selected from the list S and include 
them in R ; otherwise* make S = 0. Set X • = 1 for 
those groups of which no job is included in S. 

Set Z = Z B ±1 X , D i = { Z (s + w p ) 

_ _ J. .eS lv J. .eS 

N., 10 10 


d k -Dk * k 


J, 1 eS 

-L J 


+ S ik X i } and s k = d k - ¥ k. Recalculate 

the value of N. , ¥ J. . e R. Find max B . . N. . and 
ij 10 _ eR 10 10 

denote the corresponding job by JV-r. Set Z = Z + B-r-. N-r-. 

13 1 3 iD 

and go to Step 4. 


3.5.2 Total Profit Maximi: 
Decision Variables s 


when Machining 


>eeds are 


3. 5. 2.1 


Step 3 remains the same as given in Sec. 3. 5. 1.1 

( t i 

except the fact that instead of p.j. k » p ijk is used f° r 

/%nl /-.-i il n/v T\T T\T/^4 - rs 4“ Vs q+ -n ' ' r>cnor>Ci r» csy-> 4 - c* 4 - Vi m o vi rrn rm, 


calculating Note that 3 

production-rate unit processing time. 


represents the maximum- 
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3*5.3 Illustrative Example and Computational 


3* 5. 3.1 


Numerical example 5 presented in Sec. 3. 4. 1.5 is 
considered to demonstrate the heuristic solution procedure. 


3*5. 3* 2 Solution ; 

The following is the stepwise solution of the example 
using the heuristic approach. 

Step 1; Q 1 = 3105.00 mins, and Q £ = 2995.00 mins. 

Since Q k > d k ¥ k, go to Step 3. 

Step 3.1: Select an initial node: set n = 1, N1 = {n} = {l}. 

For node 1, set R = ^21 * ^22 9 ^23’ ^31’ ^32’ 

J n , J 42 , J43}, s = 0, Dj, = o ¥ k, s k = 240a00 min ¥ k, 

Z = 0 and X. = 1 ¥ i. 

Go to Step 3*2. 

Step 3.2: Calculate e R. Find max B ij N ij* Its 

J i0 

value is Rs. 4080.00 corresponding to job J 12* 

¥e find that ^ = 68.00 mins and 14, = 94.00 mins 
for w 12 = 10 pcs. 

Go to Step 3.3* 

Step 3.3: Since ¥ k and w-^ = l^ 5 se_t R = ^”^12^ 

= { J n» J 21* J 22’ J 23’ J 31’ J 32* J 41 ? J 42’ J 43^ * 

S = S + {J 12 } = (J 12 }» % = 0»D 1 = 68.00 mins, 

D 2 = 94.00 mins, s 1 = 2332.00 mins, and s 2 = 2306.00 mins. 
Go to Step 3.2. 
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Step 3. 


Now we recalculate the value of N. . ¥ J. . e R. It 

13 13 

is observed that max B. . N. . occurs corresponding 

Jh ..eR ^3 1 3 

to a job J 22 • We find that = 562.00 mins, 

D 2 = 473.00 mins and w £2 = 50 pcs. Go to Step 3.3. 

Since w 22 = 1 22 and \ < d^ ¥ k, set 

R = ^ J ii* J 21’ J 23* J 31’ J 32 ? J 41 s J 42* J 43^ ’ 

S = {J 12 , J 22 }» X, = 0, D- 1 = 562.00 mins, D 2 = 473.00 

mins, s^ = 1838.00 mins, and s 2 = 1927.00 mins. 


The Step 3.2 and 3.3 are followed iteratively and 
sets R & S are generated. 

R s {j^^ s J42} and 

S = {J-,,, J 


11 ? u ]_2’ u 21* u 22’ ^23 9 ^32 ? ^43^* Furthers 


, 3 " 9 33 


D 1 = 2394.00 mins, D 2 = 2386.00 mins, = 0 ¥ i and 


w 42 = 0. 


Since w 


42 


0 and D k < d v , we go to Step 3.4. 


k 


is We remove the last two jobs ^23 and 3"^ from S 
and include them in R« 


Now s set R — {^23* '^31* '^"41 , ^b2 9 ^43 ^ 9 


{J- 


J-IO. J T-Zo} » ^ 


■4 ~ 


’ll* “12* u 21* “ 22* 32 

Z = Rs. 2210.00, D x = 1504.00 mins, 

D 2 = 1536.00 mins, s x = 896.00 mins and 

s 2 = 864.00 mins. 
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Recalculate N. . ¥ J. . e R. Find max B. . N. . . Its 
1J 1J J. .eR 10 13 

value is Rs. 828.00 corresponding 1 ^ job J ^ . 

Set *Z = Z + = Rs. 3038.00 and go to Step 4. 

Now we follow B & B procedure from Step 4 onwards as 
discussed in Sec. 3. 4. 1.4. The final product mix obtained is 
as shown in Table 12. 

3 . 5 . 3 . 3 C omput a tional experience : 

It was observed that the solution of the numerical 
example 5 obtained bjr heuristic and optimizing procedures turns 
out to be the same. The heuristic and optimizing solution 
procedures required a CPU time of 0.20 secs, and 0.26 secs, 
respectively for solving this example ''nDEC-1090 computer 
system. Number of problems were solved using both the approaches. 
In each case, it was observed that the heuristic approach requires 
less computational effort as compared to the optimizing procedure. 

3.6 VARIATIONS OF MULTI-STAGE PROBLEM t 

In Sec. 3.4 and 3.5, we presented mathematical 
models and solution procedures for the multi-stage production 
systems and the objective was to maximize the total profit. The 
mathematical models and solution procedures can be suitably 
modified for maximization of total number of units produced. 

The term B. i.e., the per unit profit for job J. is 

1 J 

replaced by 1. When Ns is put equal to 1, the multi-stage 
problem becomes single stage problem. 
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The multi-stage algorithms given in Secs. 3.4- and 
3.5 were used to solve the single stage numerical examples 
lj2 s 3 and 4 and it was observed that the solutions obtained 
were same as obtained by algorithms given in Sec. 2.4. 

The solution procedure given in Sec. 3.4.2 is based 
on the fact that the unit processing times of the jobs are 
determined by machining speeds which are treated as decision 
variables. However , there may be multi-stage problems in 
which the unit processing time for some jobs are prespecified 
while for the remaining jobs the unit processing time depend 
upon the machining speeds which are also to be optimally 
selected. 

Initially, we assume that the jobs for which the 
machining speeds are not prespecified, are processed using 
minimum unit production time machining speeds. The machine 
loading and product mix is determined using the algorithm 
given in Sec. 3.4.1. Using this product mix, the time required 
on each production stage to process the selected jobs for 
which the unit processing time is not prespecified, is calcu- 
lated and subtracted from the available processing time for that 
stage. The remaining available processing time at each stage 
is used to obtain the optimal machining speeds for the various 
jobs (for which the processing times are not prespecified) to 
be processed on that stage. 
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The multi-stage problem is now converted into 
multiple single stage problems (the number of single stage 
problems considered being equal to the number of stages in 
the multi-stage problem) . The optimal machining speeds are 
determined using the basic concepts given in Sec. 2. 4. 1.2. 



CHAPTER IV 


CONCLUSIONS AND SCOPE FOR FURTHER WORK 

4.1 CONCLUSIONS : 

In this thesis, we have developed mathematical models 
and solution procedures for the machine loading and product 
mix decision problems for single and multi-stage production 
systems with limited available production times at each stage. 
The objective was to maximize either the total profit or the 
total amount of production. Two important variations of the 
problem, namely, the unit production time is prespecified and 
machining speeds are decision variables, were also considered. 

For each one of the above stated variation of the prob- 
lem, an optimizing algorithm which utilized the basic frame- 
work of branch and bound procedure, was developed. In order to 
reduce the computational effort, heuristic algorithms were 
also developed for each of the problem variations. Sample 
problems were solved using both the optimizing and heuristic 
algorithms. 

4.2 SCOPE FOR FURTHER WORK: 

In order to account for realistic production situations, 
the single and multi-stage models developed in this thesis need 
to be extended to account for the following situations: . 
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i) Each production stage has one or more machines. 

ii) Machines are capable for processing more than one unit 
at a time. 

iii) The jobs may have technological ordering for the opera- 
tions to be performed on various production stages. 

iv) The jobs may have due dates associated with them. 

In Section 2.4. 1.2, a technique was suggested for deter- 
mining the step length for decreasing or increasing the value 
of p. . Better methodology^ needs to be developed for appro- 
priate selection of step lengthj since this will considerably 
reduce the computational effort involved in the use of various 
algorithms. 
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